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Preface 

In  this  study,  the  authors  have  attempted  to  lay  the  groundwork  need¬ 
ed  to  determine  the  feasibility  of  a  trajectory-match  technique  for  inter¬ 
cepting  intercontinental  ballistic  missiles.  In  Chapter  II,  the  ballistic 
trajectory  equations  end  rocket  equations  of  motion  are  developed  in  de¬ 
tail,  since  the  authors  feel  that  many  readers,  who  may  be  interested  in 
I.Q3M  defense,  may  not  be  completely  familiar  with  these  concepts. 

We  found  the  boost-phase  analysis  in  Chapter  IV  to  be  the  most 
challenging  aspect  of  the  problem.  As  far  as  we  know,  we  have  developed 
a  new  technique  for  obtaining  good  approximations  to  the  thrust  and 
weight  of  boost-rockets  required  for  boosting  a  given  payload  onto  a 
specified  trajectory. 

Figures  are  used  liberally  throughout  the  report  so  that  the  reader 
can  more  easily  follow  the  geometry  involved  in  some  of  the  more  compli¬ 
cated  mathematical  derivations.  There  is  also  a  list  of  definitions 
included  on  page  vl  to  explain  some  terms  which  are  peculiar  to  this 
particular  study.  The  computer  programs  used  are  included  in  Appendix  A, 
so  that  other  investigators  nay  reproduce  any  of  the  results  or  expand 
this  study,  if  they  so  desire. 

Ve  wish  to  thank  Captain  H.G.  Pringle  for  suggesting  this  thesis 
topic,  and  for  his  many  helpful  ideas  throughout  the  investigation. 

The  authors  are  also  indebted  to  their  wives  and  children,  for 
their  patience,  understanding,  and  many  sacrifices  which  have  greatly 
contributed  to  the  completion  of  this  work. 

H.  A.  C. 

C.A.M. 
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Definitions 


Interceptor  Vehicle  The  final  payload  consisting  of  the 
vehicle  structure,  and  equipment  necessary  for  detection  and 
destruction  of  the  enemy  warhead. 

Velocity-Match  Rocket  An  assembly  consisting  of  the  interceptor 
vehicle,  rocket  engine,  fuel  tanks,  eto. ,  needed  to  accelerate  the 
interceptor  vehicle  to  its  final  velocity.  The  velocity  match 
rocket  is  the  payload  for  the  boost-phase  rooket. 

Interceptor  System  The  structure  consisting  of  the  boost-phase 
rocket  and  the  velocity-match  rocket. 

Excess  Time  The  difference  between  the  flight-times  of  the 
•  * 

enemy  warhead  and  the  interceptor  vehicle  from  their  launch  sites 

4 

to  the  intercept  point. 

vS*  - 

Sub-Bocket  An  assembly  consisting  of  the  payload  and  the 

.V 

^'rocket  stages  not  yet  jettisoned. 
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cross  sectional  area  of  rocket  stage 

aiming  point 

speed  of  sound 

height  above  the  earth 

coefficient  of  drag 

effective  mass  ratio  of  "boost-phase  sub-rocket 
effective  mass  ratio  of  velocity-match  rocket 
drag  force 

total  energy  per  unit  mass 
energy  ratio 

eccentricity  of  elliptical  trajectory 
functional  notation 
acceleration  of  gravity 

2 

proportionality  constant  (lb  mass  per  lb  force)(ft  per  sec  ) 

angular  momentum  per  unit  mass 

specific  impulse 

intercept  point 

kinetic  energy  per  unit  mass 

number  of  stages 

launch  site 

initial  3ub-rocket  mass  (lb^),  numerically  equal  to  sea  level 

weight 

mach  number 
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m  stage  mass  (lb  ) 

m 

nip  interceptor  vehicle  mass 

nip  interceptor  rooket  mass 

n  structural  factor 

r  length  of  radius  vector  from  center  of  the  earth 

p 

I  radius  at  intercept  point 

p 

la  radius  at  aiming  point  on  corrected  intercept  trajectory 
r^  radius  to  interceptor  trajectory 

r^  radius  to  enemy  warhead  trajectory 

S  distance 

S'  distance  travelled  by  one-stage  rocket,  starting  from  rest  in 

field-free,  airless  space 
thrust  force  of  the  jth  stage 

T*  flight  time  from  apogee  along  ballistic  trajectory 

Tx  excess  time  (See  list  of  definitions) 

Th.  ratio  of  thrust  of  jth  stage  to  weight  of  interceptor  vehicle 

J 

t  time 

u  rocket  exhaust  velocity 

V  potential  energy  per  unit  mass 

v  velocity 

range  over  earth's  surface 
x  horizontal  distance 

y  vertical  distance 

a  angle  between  thrust  vector  and  velocity  vector 

p  thrust  to  initial  weight  ratio  of  boost-phase  sub-rocket 

thrust  to  initial  weight  ratio  of  velocity  match  rocket 
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Y 

r 

6 

€ 


V 

v 


+  1 


4 

<t> 

<*> 

4> 


Lla 


mo 

moa 


P 


f 

V'u 

X 


flight  path  angle  (angle  between  velocity  vector  and  local 
horizontal) 

flight-path  angle  at  termination  of  tip-over  maneuver 

angle  between  aiming  point  and  intercept  point 

angular  difference  between  launch  site  and  burnout  point  in 

rectangular  coordinates 

gravitational  constant  of  the  earth 

marimim  acceleration  in  sea-level  wg’ s" 

burnout  angle  in  rectangular  coordinates 

angle  between  thrust  vector  and  local  horizontal 

position  angle  between  radius  vector  and  major  axis 

burnout  range;  angle  of  corrected  interceptor  trajectory 

angle  from  launch  site  to  intercept  point  for  theoretical 

intercept  trajectory 

angle  from  launch  site  to  intercept  point  for  corrected 
intercept  trajectory 

interceptor  launch  site  position  relative  to  warhead  apogee 

semi-range  angle  of  theoretical  intercept  trajectory 

semi-range  angle  of  corrected  intercept  trajectory 

atmospheric  density 

stage  mass  ratio  (boost  phase) 

overall  mass  ratio  (boost  phase) 

overall  mass  ratio  (velocity  match  rocket) 

thrust  angle  above  horizon  in  rectangular  coordinates 
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Abstraot 

In  this  thesis,  a  mid-course ,  trajectory-match  intercept  of 
intercontinental  ballistic  missiles  is  analyzed.  The  intercept 
is  accomplished  by  using  a  multi-stage  rocket  to  boost  an  inter¬ 
ceptor  vehicle  onto  a  ballistic  trajectory  in  the  same  direction 
of  flight  as  the  warhead  trajectory,  so  that  the  two  trajectories 
are  almost  tangent  at  a  selected  intercept  point.  Near  the 
intercept  point,  a  final  rocket  is  ignited  to  accelersti  the 
vehicle  and  match  the  warhead's  position  and  velocity. 

A  simplified  analysis  is  performed  for  the  case  of  a  co -planar 
intercept  over  a  spherical,  non-rotating  earth.  Intercept 
trajectories  are  determined  from  the  trajectory  geometry  and  the 
rocket  equations  of  motion.  An  approximate  method  is  derived  to 
estimate  the  thrusts  and  rocket  weights  required  to  boost  an  inter¬ 
ceptor  vehicle  onto  these  trajectories.  Equations  are  programmed 
for  an  IBM  1620  digital  computer  to  obtain  estimates  of  thrusts, 
weights,  maximum  accelerations,  and  the  difference  in  the  times  of 
f  flight  of  the  warhead  and  interceptor.  Various  combinations  of 
intercept  points  and  interceptor  launch  sites  are  used  to  determine 
the  first  stage  thrust  required  per  pound  of  interceptor  vehicle  for  in¬ 
tercept  of  low  angle,  minimum  energy,  and  high  angle,  6000  and  4000 
nautical  mile  range  warhead  trajectories.  The  effects  of  specific 
impulse,  structural  factor,  and  thrust  to  initial  weight  ratio  upon 
overall  thrust,  time,  and  acceleration  are  also  determined.  It  is 
found  that  there  is  an  optimum  value  of  thrust  to  weight  ratio  which 
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results  in  minimum  overall  thruBt  requirements. 

A  numerical  solution  to  the  equations  of  motion  is  performed  on 
the  computer  for  the  "boost  phase.  For  the  case  considered  it  is 
found  that:  the  actual  "burnout  velocity  is  within  three  percent 
less  than  desired  and  thrust  is  within  eleven  percent  of  that 
required,  thereby  proving  the  validity  of  the  original  "boost  phase 
approximations. 

The  computer  results  are  presented  graphically,  and  indicate 
that  this  intercept  technique  is  theoretically  possible  for  either 
a  manned  or  unmanned  interceptor  system.  However,  a  much  broader 
investigation  would  be  required  to  definitely  determine  its  useful¬ 
ness  for  defense  aga?^st  ICBM's. 
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AN  ANALYSIS 
OF  A 

TRAJECTORY  AND  VELOCITY  MATCH  TECHNIQUE 
FOR  INTERCEPTING 

INTERCONTINENTAL  BALLISTIC  MISSILES 

I.  Introduction 

General  Considerations 

The  high  speed  and  short  time  of  flight  of  the  intercontinental 
ballistic  missile  have  made  the  problem  of  finding  an  effective  ICBM 
defense  system  extremely  complex.  To  successfully  intercept  an  ICBM, 
an  interceptor  vehicle  must  be  positioned,  at  or  near  the  same  i'Ant  in 
space  as  the  ICBM,  at  the  same  time.  The  interceptor  must  also  be 
able  to  distinguish  the  warhead  from  decoys  and  destroy  it. 

One  system  which  is  presently  in  the  development  and  test  stage  i  - 
limited  to  intercepting  enemy  missiles  during  the  re-entry  phase  of  their 
trajectories.  For  this  type  of  intercept,  the  interceptor  vehicle 
approaches  the  warhead  praotioally  "head  on",  resulting  in  extremely  high 
relative  velocities.  This  places  very  stringent  guidance  and  accuracy 
requirements  upon  the  system.  Furthermore,  this  system  is  limited  to 
pomt-aeieaae  of  relatively  small,  strategic  target  areas.  Since  tin.  i4» 
practically  a  "last  ditch"  defense,  the  question  has  arisen:  Can  a 
successful  intercept  be  made  during  the  boost  phase  of  the  warhead 
trajectory,  while  the  rochet  motors  are  still  burning,  or  during  the 
mid-course,  free-f light  phase,  prior  to  re-entry? 
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It  would  be  difficult  to  intercept  a  missile  during  the  boost  phase 
because  of  the  very  short  warning  time  available,  and  the  usual  lack  of 
friendly  defense  bases  in  the  vicinity  of  enemy  launch  sites.  Thus, 
serious  disadvantages  also  exist  for  this  type  of  intercept. 

Now  let  us  consider  the  possibility  of  a  mid-course  intercept,  c<e£. 
a  method  of  reducing  the  relative  velocity  between  the  interceptor  and 
the  warhead.  This  can  be  accomplished  by:  launching  an  interceptor 
vehicle  on  a  ballistic  trajectory  in  the  same  direction  of  flight  as 
the  warhead;  making  this  intercept  trajectory  nearly  tangent  to  the 
warhead  trajectory;  and  accelerating  the  interceptor  vehicle  to  match 
the  warhead's  position  and  velocity  by  igniting  a  final  booster  rocket 
r.s:-.,r  the  intercept  point.  This  method  offsets  the  disadvant  gcs  1  ihc-  rt 
in  the  re-entry  and  boost  phase  intercepts.  Available  warning  time  would 
not  be  as  great  as  for  the  re-entry  intercept,  but  terminal  guidance  re¬ 
quirement  a  should  not  be  nearly  as  critical  due  to  the  lor  relative  -relo~ 

J x near  the  intercept  point.  Thus,  the  mid-course,  trajectory-match 
intercept  appears  to  be  a  promising  solution  to  the  problem  ul  iCBh  5;  ,-~.-nse. 

Statement  of  the  Problem 

The  purpose  of  this  thesis  is  to  analyze  a  mid-course,  trajectory- 
match  5.ntercept  of  an  I  CBM  warhead.  This  analysis  consist*  of  determining 
inhere  action  geometry,  the  required  thrust  and  weight  of  -V  . 
esptor  system,  and  the  difference  between  the  flight  times  of  the  ICBM 
warhead  and  the  interceptor,  from  launch  site  to  point  of  intercept. 

Since  for  this  type  of  intercept  there  is  a  low  relative  velocity 
at  the  intercept  point,  there  are  possible  advantages  to  having  a  man 
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in  the  interceptor  vehicle.  These  include  man' a  ability  to  mate  decieit 
and  evaluate  data  from  abort  range,  decoy  discrimination  equipment. 
Therefore ,  the  analysis  includes  a  consideration  of  the  maximum 
acceleration  imposed  upon  this  system,  and  sin  effort  is  made  to  keep 
acceleration  within  human  tolerances. 

Analysis  of  the  Problem 

The  analysis  of  an  ICBM  intercept  over  the  non-spherical ,  rotating 
earth,  by  an  interceptor  system  of  unknown  performance  would  be  extremely 
complex.  However,  this  analysis  can  he  greatly  simplified  by  assuming 
that  the  earth  is  spherical  and  non-rotating,  and  that  the  warhead  and 
interceptor  trajectories  are  co-planar.  Although  these  assumptions  are 
unrealistic,  the  basic  concepts  and  techniques,  when  established  fox  tluu 
-necial  case,  could  be  extended  to  the  general  intercept,  .wii  -» 
^resent  investigation  show  promise- 

Interoept  Analysis.  The  first  step  is  the  detersa?  na.  io->  'p  th'i 
geometry  of  the  intercept  for  the  tangent  trajectory.  The  equations 
are  developed  to  allow  variable  inputs  of  warhead  trajeotory,  intercept 
point  and  interceptor  launch  site. 

Next,  the  dynamics  of  the  velocity-match  phase  are  considered.  It 
is  shown  that  the  equations  of  motion  of  a  rocket  in  field  free  space  can 
ve  considered  to  be  a  perturbation  to  the  theoretical  free-1‘3  xght 
trajectory,  if  gravity  is  assumed  constant  during  the  velocity-match 
phase.  The  required  free-flight  trajectory  is  determined  by  correcting 
the  theoretical  trajectory  for  the  effect  of  this  perturbation. 

Boost  Analysis.  With  the  required  intercept  trajectories 
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established.,  the  next  problem  is  to  determine  the  booster  requirements 
to  place  the  interceptor  vehicle  on  this  trajectory.  These  require- 
ments  must  be  based  upon  the  weight  of  the  velocity-match  rocket  which 
in  turn  depends  upon  the  weight  of  the  interceptor  vehicle.  Since  this 
is  an  initial  investigation,  the  weight  of  the  interceptor  vehicle  is 
unknown.  A  vehicle  weight  could  be  assumed,  but  to  permit  greater 
flexibility,  non-dimensional  weight  ratios  are  used  to  determine  rocket 
weight  and  thrust  requirements  per  pound  of  interceptor  vehicle.  The 
boost  phase  analysis  is  further  complicated  by  the  large  number  of 
variable  rocket  parameters.  The  thrust,  weight,  specific  impulse  and 
structural  factor  of  each  stage,  as  well  as  the  number  of  stages  must  be 
established  to  continue  the  analysis.  Furthermore,  the  exact  sire  \~J 
shape  of  the  interceptor  system  and  the  velocity  profile  during  bo^t 
all  unknown.  Therefore  aerodynamic  drag  during  the  boos+  phase  oar  r~.+ 
be  exactly  predicted.  These  complications  make  a  closed-form  solv.ti  -- 
to  the  boost-phase  equations  of  motion  unfeasable. 

The  following  assumptions  are  made  to  obtain  an  approximate  solution 
to  these  equations: 

(l)  The  stiuctual  factor  and  specific  impulse  are  the  same  for  each 
stage. 

{2)  Drag  is  negligible  compared  to  thrust. 

(3)  Gravity  is  constant  in  magnitude  and  direction. 

(4)  Vertical  launch  and  tip-over  maneuver  can  be  neglected. 

(5)  An  inertial,  rectangular  coordinate  system  can  be  used. 

(6)  Thrust  for  each  stage  is  constant. 
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(7)  The  thrust  direction  is  the  seme  (constant  direction  relative 

to  inertial  spaoe)  for  each  stage. 

A  further  simplification  is  made  by  making  the  thrust  to  initial  weight 
ratio  the  same  for  each  sub-rocket. 

With  the  use  of  the  above  assumptions  and  simplifications,  a  ev.i, 
of  equations  are  obtained  which  give  the  thrust  and  weight  ratios 
required  to  achieve  a  desired  velocity  and  direction  in  spaoe.  These 
equations  are  then  combined  with  the  interoeptor  trajectory  equations 
in  a  program  designed  for  use  on  the  IBM  1620  digital  computer.  By 
means  of  this  program  the  thrusts  and  weights  of  the  booster  rockets, 
the  flight  times  from  interoeptor  launch  site  to  intercept  point,  and 
maximum  accelerations  are  obtained. 

A  numerical  integration  is  then  performed  on  the  1620  computer 
,w  check  the  validity  of  the  boost-phase  approximations.  For  inc 
•ypical  intercept  situation,  it  is  found  that  the  burnout  velocity  oi 
tue  booster  rockets  differed  from  that  desired  by  approximately  three 
per  cent. 

Effect  of  Variables.  The  large  number  of  variables  require  that 
this  investigation  be  limited  to  a  few  of  the  many  possible  enemy 
warhead  trajectories. 

Two  typical  ICBM  ranges  (6000  and  4000  nautical  miles)  are  selected 
Heaver,  this  must  be  further  limited  since  for  any  particular  varhoa,' 
range  there  are  an  infinite  number  of  possible  trajectories,  characterised 
by  the  height  of  their  apogees.  In  an  attempt  to  adequately  consider 
possible  intercept  conditions,  a  high,  medium  (minimum  energy),  and  low 
altitude  trajectory  is  considered  for  eaoh  range. 
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The  effect  of  etructual  factor,  specific  impulse,  and.  thrust  to 
weight  ratio  on  thrust  requirements,  excess  time,  and  maximum 
acceleration  is  investigated  for  a  few  typical  intercepts.  Finally, 
reasonable  values  of  specific  impulse  and  structural  factor,  and  ft- 
optimum  value  of  thrust  to  weight  ratio  are  used  to  investigate  the 
effect  of  varying  the  intercept  point  and  interceptor  launch  site  on 
thrust  and  excess  time. 
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II.  Basic  Theory 

The  Free-Flight  Missile  Trajectory 

Basic  Equations.  The  motion  of  a  Ballistic  missile  in  a  vacuut. 
under  the  influence  of  a  central  force  field  such  as  the  earth*  0  yravi-- 
tational  field,  can  he  analyzed  hy  considering  the  missile  to  he  a 
point  mass  in  a  polar,  inertial  coordinate  system  (Ref  13:2-20). 


Let  i  he  a  unit  vector  in  the  r  direction  and  j  a  unit  vector  per¬ 
pendicular  to  r  as  shown  in  Pig.  2-1.  The  derivatives  of  the  unit 
vectors  with  respect  to  4>  are 


dt 
d  t 


and 


Then  the  derivatives  with  respect  to  time  are 


di  _  di  dt 
dt  d4>  ait 


and 


—  -  d>T 

dt  dt  dt 
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Now  the  position  vector  of  the  missile  is 


Its  velocity  is 


r=  re 


r  =  r  i  +  r 


r  -  ri  +  r<P  j 


and  its  acceleration  is 


or  r  =( r  -  r<pZj  L  -h{2.rcp  -h  rep)  j 

How  ^r(r2<p)  ~  2rr<p  -h  rzcp  =  r(Z  r<j>  +  r<p) 


~r~  Jt(r2<p)  =  2rj>  +  rp> 


(2-1 } 


(2-2) 


(2-3) 


Substituting  Eq  (2-3)  into  Eq  (2-2), 


r-(r-ri‘)i  +[yr  3r(r‘i)] j 


Now,  by  applying  Newton's  Second  Law  of  Motion,  P  ■  ma 


where 


£7  m  — 

F=-~pr-  l 


m  is  the  mass  of  the  missile 


p  is  the  earth's  gravitational  constant 


we  get 


Mm  ~ 


L~  m[(r-r<pz)i  -h(yr  ^t[rz<p]), 'jj 


Then  by  equating  the  j  components  of  Eq  (2-5) 

I  d  /  z  ;  \  - 


~’dT(r^*)  =  ° 


(2-5) 


Now,  since  r  is  not  infinite. 
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lit  (r&4>)=  0 


or  racfi>  —  "a  constant " 

But  r  <b  is  the  polar  coordinate  expression  for  the  angular  momentum  ; 
unit  mass.  Therefore,  the  angular  momentum  of  the  missile  is  a  constant 
during  the  free-f light  phase.  Let  us  call  the  angular  momentum  par  ~-r  1 
mass,  h.  Then 

h=rz4>  (2-6) 

Now,  equating  the  I  components  of  Eq  (2-5) > 

-  ~  r  -  r<j>*  (2-7) 


In  order  to  solve  Eq  (2-7),  let  us  define  a  dummy  variable, 


Then,  in  terms  of  p,  Eq  (2-6)  becomes 


h=r*4‘  =  jz 


or 


Also 


4>  =  ph 


But 


d<f> 

_ L 

d4>\  P  J  ~  P*  d<f> 

Combining  Eqs  (2-9)  and  (2-10), 


Then 


d<p 


(2-B) 


(2-9) 


(2-10) 


/  -  \ 
V  ^  *  ' 
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(2—1?) 


or 


dZp 

d<t>9 


Substituting  Eqs  (2-8  and  (2-12)  into  Eq  (2-7),  with  r  ■  ^  ,  we  get 

2  p 

Dividing  Eq  (2-13)  by  -pa  gives 

(~M/h*)  =  -t  p  (2-14, 

Equation  (2-14)  is  a  non-homogeneous,  second  order,  linear  differential 
equation  which  can  be  solved  by  the  method  of  undetermined  coefficients 
(Ref.  14*36-41). 

The  solution  of  this  equation  is 


p  =  jp:  -h  A  COS  (4>~  <t>')  (2-15) 

where  A  and  <b  '  416  constants  which  must  be  determined  from  the  boundc. 
conditions. 

The  evaluation  of  these  constants  depends  upon  the  feet  that  thr 
.otal  energy  of  the  missile  must  be  a  constant  for  free-flight  in  a 
central  force  field.  Thus, 

E  -  K  +  V  -  "a  constant" 

where:  E  »  total  energy  of  the  missile  per  unit  mass 

K  -  kinetio  energy  of  the  missile  per  unit  mass 
V  -  potential  energy  of  the  missile  per  unit  mass 
By  defining  potential  energy,  V,  as  zero  when  r  is  infinite  we  can  v*.y 


that 


V=  - 


r 


Also,  kinetio  energy  is  defined  as 


or 


K=2~(rz-t-  r*<pz) 


(2-16) 

(2-17) 
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Thus 


s')  M 


E-±(t-'  +  rV)- 


Equation  (2-18)  can  be  rewritten  in  terms  of  p  by  the  use  of  Eqs  (2-3' 


and  (2-8). 


£  ~  ~k[hz(ijf)  +-  (~p)(pzhf]  ~JJ  P 


h *  rsdp\i 


2  Hd 


V  +  P  j  ~^P 


(2-...V) 


Differentiating  Eq  (2-15)  gives 


^if  =  -  A  Sin  (4>-cP') 


(2-20) 


By  using  Eqs  (2-15)  and  (2-20),  p  is  eliminated  from  Eq  (2-19»  giving 
an  expression  for  E  in  terms  of  A  and  h. 


h2Aa 

£ =  a  2hz 

Rearranging  Eq  (2-19a)  and  solving  for  A, 


Substituting  Eq  (2-2l)  into  Eq  (2-15) » 


(2-19a) 


+  *  -V  ,,  2Ehz 


h‘~  h 


p~j?['  cosc*-^‘)] 


(2-22) 


J_=  h%U _ 

P  cos(^~^) 


(2- v' 


Now  if  we  choose  our  coordinate  system  so  that  r  is  a  maximum  when 
<p  •  0,  then  4>/  m  0  and  Eq  (2-23)  becomes 
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( 2-24  J 


h*/M 

I  -  V I  +  -^fr  cos  4 


This  equation  is  similar  to  the  equation  of  an  ellipse  in 
polar  coordinates,  which  is 


t-ecosf  (2-25) 

where  1  o  the  semi-lattice  rectum  of  the  ellipse 

e  -  the  eccentricity  of  the  ellipse 
Therefore  we  can  state  that  the  equation  of  the  trajectory  of  a  ballis¬ 
tic  missile  in  free-flight  over  a  spherical  non-rotating  earth  is  an 
ellipse  with  the  center  of  the  earth  at  one  focus.  The  eccentricity  of 
the  ellipse  is  defined  by  the  equation 


and  the  semi-lattice  rectum  by 


(2-27) 


Thus  the  trajectory  of  a  ballistic  missile  is  completely  defined  if 
the  total  energy  and  angular  momentum  of  the  missile  are  known,  and  is 
given  by  the  equation 

hVM 

/-s  cos  4> 

Additional  Trajectory  Parameters.  Figure  2-2  shows  the  elliptical 
trajectory  plus  several  other  parameters  which  are  useful  for  analyzing 
ballistic  trajectories.  These  parameters  are  defined  as  follows: 

Y  -  the  angle  between  the  flight  path  and  the  local  horizontal 

v  -  theoretical  launch  angle  at  the  earth's  surface 

1  0 
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vq  -  theoretical  launch  velocity  at  the  earth's  surface 

X  -  range  over  the  earth's  surface  from  theoretical  launch 

© 

site  to  theoretical  impact  point 

r  -  radius  of  the  earth 
o 

4>  -  semi-range  angle  of  the  trajectory 

T  0 

4>  -  angle  from  major  axis  of  the  ellipse  to  the  radius 

vector  (measured  positive  clockwise) 

The  angle  of  the  tangent  to  the  trajectory,  can  k©  determined 
from  the  ellipse  equation  and  Fig.  2-3. 
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tar)  r  -  c/r  =  1  dr 
ion  r  -  -  r  d<t> 

(2-2;? 

Differentiating  Eq  (2-25), 

dr  - 1  e  sin  4> 

d<t>  (l-e  cos  <£)* 

Then 

tan  r-  ,'eW, 

1  ~  <2  COS  <p 

(2-31) 

Prom  Eq  (2-31),  the  theoretical  launch  angle  is  defined  "by  the  equation 


Tiae  of  Flight.  One  very  important  parameter  of  the  trajectory 
remains  to  he  determined.  This  is  the  time  of  flight  during  the  free- 
f light  trajectory  (Ref.  13» 17-19)*  An  inspection  of  Pig.  2-3  shows  that 
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the  incremental  area  between  two  radial  lines,  d<£  apart,  is 


dA  =  ~^( rd4>)(r)  «  rzd4> 


Equation  (2-36)  can  he  written  as 


^  .  I  id*  ,  c  . 
d  t  ~  2.  r  dt  ~  2  r  ? 


I 


hut 


=  /? 


(2-3.6) 


(2-6) 


and  so 


dA 
dt  ~  2 


=  ~5~  h 


(constant) 


(2-37) 


For  one  complete  orhit  or  revolution, 

A=~^hT  (2-38) 

where  A  •  the  total  area  of  the  ellipse 

T  -  the  period  of  revolution 

Since  the  time  rate  of  change  of  the  area  of  a  sector  of  the  ellipse  it» 
a  constant, 

T*_  A* 

T  A  (2-39) 


idie re  T*  is  the  time  required  to  travel  along  the  arc  of  a  sector  of 
area  A*. 


Then 


T  A  It 


But,  hy  rearranging  Eq  (2-38), 


A  * 

T  _ 
A 


4  r  ,  i 


2' 


Thus 


T*=fA* 


(2-41) 


How,  since  A*  is  the  area  of  any  sector  of  the  trajectory,  let  us 
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evaluate  A*  for  a  sector  from  the  apogee  radial  line  (where 
to  the  general  angle ,  4>  ,  downrange  from  apogee. 

From  Eq  (2-36)  . 

dA  * -j-  rzd<f> 

Squaring  Eq  (2-25)  gives 

Za  _ 

r  ~  (/-  e  cos  <f)*  ^A(l~e  cos  <p)z 
Combining  Eqs  (2-36)  and  (2-42) 


Then 


dA  = 


h4d4> 

2^u2(l-  e  cos  <£)* 


A>  h*  f*  d* 

Jo0-e  ccs  4>)z 


-  o) 


(2-42) 


(2-43) 


(2-44) 


This  integral  may  be  evaluated  by  the  use  of  a  table  of  integrals 
(Ref.  8t 41)42).  The  result  is 


/  h4  )f  I  \f~  e  sin4> 
2yuz/\l- e Vjf  /-e  co 5 


,  ~i Ml-*'  tan 
ian  ( — — 


Substituting  Eq  (2-45)  into  Eq  (2-41)  we  obtain  the  equation  for  the 
time  of  flight  from  apogee  to  any  angle ,  »  downrange  from  apogee. 


T 


{Mdpfc 


e  st n  4> 


e  cos  4 


-c*  tan 
/ -e 


(2-46) 


Since  the  trajectory  is  symmetric  about  the  apogee  radius  (or  nol 
axis))  Eq  (2-46)  can  also  be  used  to  determine  the  time  to  apogee  from 
any  angle  up-range  by  substituting  the  absolute  value  of  the  position 
angle  in  the  formula. 

The  time  of  flight  between  any  two  points  on  the  trajectory  is 
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X*-  T*M  -  TU) 

(2-47) 

where  <£,>  <f>^  and  the  correct  signs  are  used  for  the  position  angle , 
Minimum  Energy  Trajectory.  (Ref.  13*1924)  Previously?  we  have 
shown  that  the  trajectory  of  a  ballistic  missile  is  completely  defined 
by  the  total  energy,  E,  and  the  angular  momentum,  h.  Now  we  will 
investigate  the  effect  of  each  of  these  parameters  on  the  trajectory. 

First,  let  us  note  from  Eq  (2—48)  that  sb  the  total  energy  increases, 
the  theoretical  launch  velocity,  v,  increases. 


-L  VZ  — 
a  v° 


re 


(2-48) 


But 


h=  r0v0  cos  r0 


(2-33) 


Thus  for  a  constant  launch  angle,  the  angular  momentum  also  increases 
with  an  increase  in  total  energy. 

Rearranging  Eq  (2-28)  gives  an  equation  which  relates  the  angular 
momentum  and  eccentricity  to  the  semi-range  angle 


COS  <t>0 


Substituting  Eq  (2-26)  into  Eq  (2-49)  gives 


(2-49) 


(2-50) 


Inspection  of  Eq  (2-50)  ehows  that  as  E  increases  with  h  held 
constant,  cos  <fi  decreases  or  <fi  Q,  and  therefore  range,  increases 
(See  Fig.  2-4a). 

Examination  of  Eqs  (2-33,  48,  and  50) shows  that  for  a  constant 
launch  angle,  an  increase  in  the  total  energy  increases  both  the  range 
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and  the  angular  momentum  (Pig.  2-4h ) . 

Now  let  us  determine  the  effect  of  angular  momentum  on  the  total 
energy  required  for  a  fixed  range.  Equation  (2-50)  can  he  written  ?? 


If  there  is  a  minimum  energy  for  a  given  range,  we  can  determine  the 
value  that  h  must  have  for  minimum  energy  hy  differentiating  Eq  (2-5l) 
with  respect  to  h,  equating  the  result  to  zero,  and  solving  for  h. 
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Finally,  for  a  given  range  there  is  an  almost  unllmit-.,i 
trajectories  possible  depending  on  the  values  of  E  and  h.  For  any 
particular  value  of  E  which  is  greater  than  E^^,  there  are  two  possible 
trajectories  for  the  same  range.  These  two  trajectories  converge  to  t>? 
minimum  energy  trajectory  as  the  total  energy  approaches  E  ^  (See  Fig. 
2-4d). 

Now  that  we  have  established  the  existence  of  a  minimum  energy 
trajectory,  let  us  derive  a  few  minimum  energy  relationships.  Substi¬ 
tution  of  Eqs  (2-54)  and  (2-55)  into  Eq  (2-26)  gives  the  eccentricity 
of  the  minimum  energy  trajectory  as 

Cos  <£,_ 

(2-56) 


&  — 
' —  M  P" 


ME-  I  -h  sin  4> 

The  optimum  launch  angle  is  determined  by  substitution  of  Eq  (2-5 6) 
-to  Eq  (2-32). 


tan  x 


_  cos  4>o 


5m.£  / -f  sin  4* 

1  o 

Notice  for  the  minimum  energy  case, 


(2-5?) 


e.  tan  XMJE 


Energy  Ratio  Concept.  The  energy  ratio  (Ref.  6:  134)  is  a 
measure  of  the  energy  available  for  the  trajectory,  and  i3  defined  as, 

ra  VB 


ER  = 


If  the  range,  or  <p>  and  the  energy  ratio  are  specified,  the  launch 
angle  needed  to  achieve  this  range  is  determined.  The  equation  is 

tan  rs  =  r/rlfx  (2-58) 
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This  equation  is  used  to  specify  the  high  and  low  angle  trejec:t'  :’ic  ;■ 
which  are  investigated  in  the  latter  portion  of  this  report* 


Rocket  Equations  of  Motion 

Velocity  Equations.  The  motion  of  a  rocket  can  "best  be  describe! 
by  considering  the  case  of  a  rocket  in  field-free  space  (no  external 
forces).  Let  us  assume  that  we  have  a  rocket  in  this  environment , 
traveling  horizontally  with  velocity,  v,  and- mass,  M,  at  time,  t. 

(Fig.  2-5). 

E  >  ~*v 

Rocket  at  time *  t 

u-v  M-mt  ^  — ►  v+dv 

Rocket  at  time,  t  +  dt* 

Fi$.  2-JT 

Diagram  -for  Rocket  Equation  Derivations 


Let  t. 

o 

M 


s 

0 


% 

m 


u 

Since  no 
momentum  must 


=  the  total  burning  time  of  the  rocket 

■  the  mass  of  the  rocket  before  ignition 

=  the  mass  of  the  rocket  at  burnout 

■  the  velocity  of  the  rocket  when  t  *  0 

=  the  velocity  of  the  rocket  when  t  *  t^ 

=  the  displacement  (position)  of  the  rocket  at  x  *»  ■.’ 

a  the  displacement  of  the  rocket  at  t  =  t^ 

*  mass  rate  of  flow  (burning  rate)  of  fuel,  a  constant 

*  the  exhaust  gas  velocity  relative  to  the  rocket 
external  forces  are  acting  on  the  rocket  of  Fig.  2-5,  the 
be  a  constant  for  the  rocket  system  (including  exhaust 
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(2-69) 
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s-s.i-v.t  +  ujjnpi^dt 

(2-?^; 

u  M  T  Mm 

Hence  ,S-  so~h  V0t  +■  ut  .  In  —rr 

O  o  m  A7 

(2-75) 

For  the  burnout  conditions,  Eq  (2-75)  becomes 

Sb  -  Sa  +  V0tk  -b  S' 

(2-76) 

wh9r0  ,  ,  uMf  7  M. 

(2  7  *' 

or  5=  u[tt  --%lnD] 

(2-78) 

Equations  (2-70,  76,  and  73)  are  the  "basic  rocket  equations  which 
used  in  the  analysis  of  the  problem.  They  can  be  modified  as  necessary 
to  account  for  the  effects  of  gravity  and  atmospheric  drag. 
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Specific  Impulse  of  Rocket  Fuels.  The  velocity  and  distance 
equations  of  rocket  motion  are  both  functions  of  the  velocity  of  the  roo 
exhaust  gases.  Usually,  the  exhaust  velocity,  which  is  a  characteristi 
of  a  given  fuel,  is  not  stated  explicitly.  Instead,  the  specific 
impulse  of  the  fuel  is  used  to  define  the  exhaust  velocity  by  the 
equation  (Ref.  10  s  20) 

u  »  Is  gc  (2-79) 

where  I  ,  the  specific  impulse,  is  the  pounds  of  thrust  per  pound  of 
s 

mass  flow  per  second,  and  g  is  a  proportionality  constant  numerically 

c 

equal  to  the  sear-level  acceleration  of  gravity. 

The  derivations  which  follow  in  Chapters  III  and  IV  have 
exhaust  velocity,  u,  as  a  variable.  However,  the  computer  progrsssp  up* 
to  solve  the  equations  have  I#,  rather  than  u,  as  an  input  variable. 
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III.  Determination  of  Intercept  Trajectory 


Theoretical  Intercept  Geometry 

As  a  first  step  in  determining  the  interceptor  trajectory  let  us 
compare  two  possible  cases  as  shown  in  Pigs.  3-1  and  3-2.  In  the 
"apogee"  intercept,  the  point  of  intersection  between  the  warhead 
and  intercept  trajectories  occurs  when  the  interceptor  is  at  the 
apogee  of  its  trajectory.  On  the  other  hand,  for  the  "tangent"  inter¬ 
cept,  the  intersection  is  at  the  common  tangent  to  both  the  warhead 
and  intercept  trajectories,  and  this  point  is  "down  range"  from  the 
interceptor's  apogee.  By  comparing  Pigs.  3-1  and  3-2  it  is  seen  that 
for  the  same  intercept  point  and  interceptor  launch  site,  the  tangent 
intercept  will  result  in  a  smaller  vector  velocity  difference  between 
the  warhead  said  the  interceptor.  This  is  an  important  fact,  since 
the  weight  of  the  velocity-match  rocket  increases  exponentially 
with  the  magnitude  of  this  velocity  difference.  Therefore  this 
investigation  is  limited  to  the  tangent  intercept  technique. 

The  most  direct  approach  would  be  to  fix  the  desired  intercept 
point  and  interceptor  launch  site,  and  then  determine  the  intercept 
trajectory  parameters.  However,  this  leads  to  several  transcendental 
equations  which  can  not  be  solved  explicitly  for  the  desired  quantities. 
As  an  alternative,  the  following  indirect  approach  is  considered. 

Prom  Chapter  II,  the  equation  of  the  warhead  trajectory,  over  a 
spherical,  non-rotating  earth  is 


r"  ^ (!-e 


hi 


W 


COS 


4J) 


(3-D 
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where 


is  the  dlstanoe  from  the  oenter  of  the  earth 


h^  is  the  angular  momentum  per  unit  mass 
e^  is  the  eccentricity  of  the  ellipse 

is  the  angle  measured  from  the  warhead  polar  axis 
Similarly,  the  equation  of  the  interceptor  trajectory,  ooplanar  with 
the  warhead  trajeotory,  is 

hi _ 

^  A  *  \ 

(3-2) 


r  — 


^0-em  ccs  4>m) 


Now  at  the  point  of  intercept ,  the  radii  of  "both  the  intercept  and 
warhead  trajectories  must  he  equal ,  or 


<NX 


=  ri 


(3-3) 


where  r^  is  defined  as  the  distance  from  the  center  of  the  earth  to 
the  intercept  point.  Furthermore,  if  we  specify  that  the  two  tra¬ 
jectories  are  tangent  at  the  intercept  point,  the  path  angles  must  be 
the  same  for  both,  or 


rmI  « 


- 


(3-4) 


where  fj  ts  the  angle  measured  from  the  local  horizontal  to  the  tangent 
to  the  trajectories  (positive  counterclockwise).  This  angle,  fj*  can 
be  found  by  using  Eq  (2-31).  Therefore 


or 


Iran 


rz  = 


-e.  sin  <t>„x 
cos  4>w:z 


tan 


Xl  - 


-  fw  sin  4>~x 

l-em  cos  4mJ. 


(3-5) 


(3-5a) 


Rewriting  Eq  (3-2)  for  the  intercept  point,  and  by  using  Eq  (3-3) 


we  get 


rx 


_ h 


2 

m _ 

cos 


(3-6) 
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We  can  also  write  Eq  (3-2)  for  the  initial  conditions,  or  launoh 
site  on  the  earth's  surface,  for  which  oase 

r_ _ hi 

°  t~em  cos  <fime)  (3-7) 

where  is  the  semi-range  angle  of  the  intercept  trajeotory. 

For  a  fixed  warhead  trajectory  and  intercept  point  we  know  e^,  h^ 
<£WI»  Tj  and  Yj*  Equations  (3-5*,  6  and  7)  oontain  four  variables,  e^ 
h  ,  d>  ^  and  ep  .  Therefore,  if  we  fix  4>  ,  we  have  three 

equations  and  three  unknowns  and  the  theoretical  intercept  trajeotory. 
parameters  are  completely  determined.  To  find  these  parameters,  we 
first  divide  Eq  (3-6)  by  Eq  (3-7)  to  give 

rx  __  /  -  em  cos  4>mo 

r°  /_e«  cos  4>mx  (3_8) 

Now  solving  Eq  (3-5a)  for  e  ,  we  get 

m 

e  »  -tan  Xr _ _ 

m  SLn  ^rnx  -  tan  Tz  cos  4>mX  (2,-9) 

Substituting  Eq  (3-9)  into  Eq  (3-8)  and  simplifying, 


sin  4>mT  -  tan  7TX  cos  4>mx  +  tan  4mX  cos  4>mo 

SLn 


(3-10) 


By  rearranging  Eq  (3-10)  we  can  get 


sin  4>mx 
tan  rx 


(-*  ' ')  +  005 


=  003  4>mc 


(3-11) 


Next,  let  us  define  a  new  quantity  which  is  a  constant  for  the 
intercept  point.  Let 


/ 

tan  rz 


(3-12) 
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Substitute  Eq  (3-12)  into  Eq  (3-ll)  and  rearrange. 

cos  4>mx  +  'S  sin  4>mx  =  cos  4>mo  (3-13) 


Equation  (3-13)  contains  only  one  unknown,  <f>  ,  and  this  can  be 

solved  as  follows?  Construct  a  right  triangle  as  shown  in  Pig.  3-3. 
By  referring  to  Pig.  3-3  it  can  be  seen  that 


/  =  sin  A  (3-14) 

cos  A  (3-15) 


If  we  substitute  Eqs  (3-14)  and  (3-15)  into  Eq  (3-13)  the  result  is 


'fi+Yr sin  A  COS  4>mx  +  yV  y*1  cos  A  sin  4>mx  =  cos  +mo 


(3-17) 


which  can  be  rewritten  as 


By  solving  Eq 


sin(A  +  = 

(3-18)  for  <j>mx  , 


cos  4>„b 

V7*  (3.18) 

and  using  Eq  (3-16),  the  final  result  is 
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^mx 


sen 


(3-19) 


With  4>  now  determined,  we  can  solve  for  the  intercept  trajectory 
parameters  em,  and  h^,  hy  using  Eqs  (3-9)  and  (3-7)»  We  have  also 
indirectly  determined  the  theoretical  launch  site  of  the  interceptor 
system  since  the  angle  from  launch  site  to  intercept  point,  *£^1’  is 
just  (See  Pig.  3-2) 


Az  ~  +-  I 


(3-20) 


Velocity  Difference  Between  Warhead  and  Interceptor 

With  the  theoretical  intercept  trajectory  now  completely  defined, 
we  wish  to  find  the  difference  in  velocities  between  the  warhead  and 
interceptor,  so  that  we  may  determine  how  large  a  rocket  is  required 
to  perform  the  velocity-match,  i.e.  reduce  the  velocity  difference  to 
zero.  Ignition  of  the  velocity-match  rocket  causes  the  actual 
trajectory  to  depart  from  the  theoretical  trajectory,  and  an  approximate 
method  of  correcting  for  this  departure  is  discussed  later.  By  referring 
to  Eq  (2-33)  it  is  seen  that  the  angular  momentum  per  unit  mass  of  the 
warhead  or  the  interceptor  can  be  written  as 

Av  «  rx  Kvr  css  rx  (3‘21) 

and  / 

hm  =  rt  \zmx  cos  rz  (3_22) 

If  we  now  divide  Eq  (3-21)  by  Eq  (3-22)  we  get 

fom  _  Vmx 

h,  ~  v„  (3-23) 
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which  can  he  solved  for  v  to  give 

mi 


y  —  \J  _l£! 

**  wr  h„ 


The  velocity  difference,  A  v  is  given  hy 


^  V  =  ~ 


(3-24) 

(3-25) 


By  making  use  of  the  fact  that  v  „  and  v  are  co-linear,  and 

nx  wx 

substituting  Eq  (3-24)  into  Eq  (3-25)  we  find  the  magnitude  of  the 
velocity  difference  to  be 

Av/=  v„0-t) 

It  is  interesting  to  note  that  if  the  angular  momentum  of  the 
interceptor  approaches  that  of  the  warhead,  A  v  approaches  aero,  but 
with  the  intercept  conditions  we  have  imposed,  this  occurs  only  when  the 
two  trajectories  are  identical  ...  a  highly  impractical  oase. 


Correction  to  Theoretical  Intercept  Trajectory. 

In  order  to  correct  the  theoretical  intercept  trajectory  for 
the  effect  of  the  velocity-match  rooket,  we  must  make  use  of  the  rocket 
equations  of  motion  which  were  developed  in  Chapter  II.  Equation 
(2-70)  can  be  written  in  vector  notation  as 

7b  =  V„  -h  uln  D  (3-27) 

Equation  (2-76)  can  also  be  written  in  veotor  notation  as 

-  ■?.  +  v.  %  +  S'  (3-28) 

where  s'  is  now  defined  by 


(3-29) 
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If  we  let  our  coordinate  system  coincide  with  the  point  where  the 
velocity-match  rocket  is  ignited,  then  bq  »  0  and  Eq  (3-28)  becomes 

=  Vctb+S'  (3-30) 

Now  if  we  re-write  Eqs  (3-27)  and  (3-30)  and  assume  that  the  motion 
takes  place  in  a  gravitational  field,  where  g  is  constant  in  magnitude 
and  direction,  these  equations  are  modified  to 

Vk=V0-buln  D  (3-31) 


and 


5.  =  V.tk  ~hs'  +T5t 


(3-32) 


Let  us  now  refer  to  Fig.  3-4  and  assume  that  our  interceptor 
rocket  is  in  free-flight  on  the  theoretical  interceptor  trajectory. 

If  we  did  not  fire  the  velocity-match  rocket,  we  would  coast  to  the 
intercept  point,  IP,  and  arrive  there  with  velocity  equal  to  vni  . 

But  we  want  to  match  the  warhead's  velocity,  v^r  ,  at  the  intercept 
point.  Therefore  the  velooity-oatch  rocket  must  he  ignited  some 
time  before  we  get  to  point  IP.  Let  us  assume  that  the  velocity-match 
rocket  has  a  burning  time  equal  to  the  free-flight  time  from  point  A 
to  IP.  Then,  when  the  rocket  arrives  at  point  A  with  velocity  vq, 
the  velocity-match  rocket  is  fired  so  that  the  exhaust  velocity  vector, 
u,  is  parallel  to  the  trajeotory-tangent  at  point  IP.  If  the  effective 
mass  ratio,  D,  of  the  velocity-match  rocket  is  of  the  right  magnitude, 
then  the  burnout  velocity,  v^ ,  will  be  approximately  equal  to  v  ,  as 
seen  in  Fig.  3-4.  However,  Fig.  3-4  also  shows  that  the  position  at 
burnout  is  point  B  rather  than  that  desired  (IP). 

Suppose  we  correct  this  discrepancy  by  the  method  shown  in 
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Fig.  3-4 

Vector  Diagram  of  Velocities  and  Displacement 
■for  Theoretical  Intercept  Trajectory 
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Fig*  3-5*  First  we  extend  the  tangent  line  ’’up-range*'  to  an 
aiming  point,  AP,  so  that  AP-IP  is  equal  to  IP-B.  Next  we  adjust 
the  trajectory,  so  that  a  free-flight  interceptor  will  arrive 
tangent  to  point  AP  with  the  same  velocity,  vffiX  ,  that  existed  at 
point  IP  of  the  theoretical  trajectory.  Now,  if  the  velocity-match 
rocket  is  fired  at  some  point  A’ ,  where  A*  is  chosen  so  that  the 
free-flight  time  from  A'  to  AP  is  equal  to  t^,  the  burnout  position 
will  be  at  IP.  The  actual  boosted  trajectory  will  be  along  the 
dotted  line  in  Pig.  3-5* 

If  we  assume  that  the  burning  time  is  short,  we  can  say  that 

ffl  «lrx\  (3-33) 


Then  the  change  in  direction  and  magnitude  of  gravity  between  points 
A'  and  IP  will  be  extremely  small,  and  by  referring  to  the  velocity 
vector  diagram  in  Pig.  3-5  we  can  say  that 


Vo  +  9  **  Kwx 


(3-34) 


uln  D  -  v \  -  gtk  -  vc 
uln  D=v,-Vmz 


(3-35) 

(3-36) 


However,  we  want  v^  to  be  equal  to  v  ,  and  since  u,  and 

v  are  all  parallel  we  can  write  Eq  (3-36)  as 
mi 


uln  D  =  vwz  -  vmI  «  A  v 


(3-37) 


where  4v  is  given  by  Eq  (3-26).  By  solving  Eq  (3-37)  for  D  we 


find  that 


D  =  e 


Av/u 


(3-38) 
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Fi  cjf.  3-5 

Vector  Diagram  oi  Velocities  and  Displacement 
•for  Corrected  Intercept  Trajectory 
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In  Chapter  IV  it  will  he  shown  that  if  we  speoify  the  thrust  to 
initial  weight  ratio,  0,  of  a  rocket,  we  can  find  s'  such  that 

s'=  H u,  P,  D) 

With  D  and  s'  known  we  can  calculate  the  corrected  intercept  trajectory 
parameters. 

By  referring  to  Fig.  3-6  it  is  seen  that 

4-  F=90°-'TX  (3-39) 

Note  that  Yj  is  negative  for  intercept  points  down-range  from  the 
warhead  apogee.  We  can  also  show  that 


tan  6  — 


s'sin  F 

nT  -  s'cos  P 


(3-40) 


or  using  Eq.  (3-39) 


,  r  >5  COS 

tan  o  -  — — -r~. — — 
rx  -  s  sin  irx 


(3-41) 


Then  with  6  known  we  can  show  that  the  flight-path  angle  at 


point  AP  is 

-Ji -s 


(3-42) 


where  6  is  always  a  positive  angle  and  the  correct  signs  must 
he  used  for  yja  and  yj*  By  using  the  law  of  sines  we  find  that 


r  -  £lf  LH  H 
■r‘,  sin  S 


(3-43) 


or 


_  _  s'  cos  Ik 

~  sin  6 


(3-44) 


With  Yt  i  rT  »  and  v  _  known  at  the  point  AP  we  can  prooede  to 
'la  la  ml 

find  the  trajectory  parameters. 

Re-writing  Eq  (3-22)  for  the  corrected  trajectory  at  the  point 
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The  angle  from  the  corrected,  interceptor  launch  site  to  the  aotual 
intercept  point  is 

=  +  ftnxa  +  &  (3-51) 

where  6  is  given  hy  Eq  ( 3—41 ) . 
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In  this  chapter,  first  we  derived  the  theoretical  intercept 
trajectory  equations  hy  fixing  the  warhead  trajectory  and  the 
intercept  point,  and  imposing  the  condition  that  at  the  intercept 
point  the  warhead  and  intercept  trajectories  had  the  same  radius, 
r^,  and  the  same  path-angle,  y^.  By  assuming  a  short  burning 
time  for  the  velocity-match  rocket,  and  that  g  was  constant  in 
the  vicinity  of  the  intercept  point,  we  then  corrected  the 
theoretical  intercept  trajectory  for  the  effect  of  the  rocket  dynamics, 
and  obtained  the  parameters  for  the  corrected  intercept  trajectory. 

We  also  determined  the  effective  mass  ratio  required  for  the 
velocity-match  rocket.  With  this  trajectory  completely  determined, 
and  the  effective  mass  ratio  of  the  velocity-match  rocket  known,  we 
can  now  consider  the  boost-phase  of  the  intercept. 
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IV.  “Boost-Phase  Analysis 

In  Chapter  III  it  was  shown  that  if  the  warhead  trajectory, 
intercept  point,  and  semi -range  angle  for  the  interceptor  were 
fixed,  a  unique,  free-flight  trajectory  is  determined.  We  also 
found  the  required  mass  ratio  for  the  velocity-match  rocket.  In 
effect,  we  are  now  at  the  point  where  we  have  a  "payload*' ,  and  a 
desired  trajectory,  and  must  find  the  "booster  rocket  requirements 
needed  to  get  this  payload  onto  this  trajectory. 

There  are  an  infinite  number  of  points  where  we  can  enter  this 
intercept  trajectory,  and  for  each  point  the  required  burnout 
velocity  and  direction  is  different.  This  point  must  also  be  at 
some  altitude  above  the  earth  because  of  the  time  required  to 
accelerate  to  the  final  desired  boost-phase  burnout  velocity.  If 
we  also  consider  the  fact  that  most  rocket  launchings  consist  of  a  vertical 
flight  segment,  a  segment  at  constant  angle  of  attack,  a  segment  at 
zero  lift,  and  a  constant  attitude  segment,  it  can  be  seen  that  there 
are  many  possible  flight  paths  from  the  launch  site  to  the  burnout 
point.  The  problem  is  further  complicated  by  atmospheric  drag, 
variable  gravity,  and  the  variation  of  rocket  thrust  with  altitude. 

Thus  there  does  not  appear  to  be  a  neat,  closed-form  solution  to  the 
problem,  although  many  approximate  methods  have  been  devised  (Ref  12). 

In  WADC  TR  57-724  (Ref  2)  a  method  for  the  preliminary 
design  and  optimization  of  a  long-range,  ballistic  missile  is 
outlined.  Basically,  it  is  a  trial  and  error  method,  and  while  it 
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is  well  suited  to  the  design  of  one  particular  missile,  it  does 
not  lend  itself  to  this  investigation,  shore  we  have  many 
trajectories  to  consider. 

Appro ilmat ions 

To  proceed  with  the  problem,  let  us  make  some  assumptions 
which  enable  us  to  develop  a  method  for  approximating  the  required 
thrusts  and  weights,  and  the  available  warning  time  for  the 
interceptor  system.  First  we  assume  that  the  boost-phaee  is  a  small 
portion  of  the  intercept  trajectory,  allowing  us  to  consider  gravity 
as  a  constant  vector.  This  enables  us  to  choose  an  inertial, 
rectangular  coordinate  system  on  the  earth’s  surface  near  the 
launch  site,  as  in  Fig.  4-1*  Next  we  neglect  the  vertical,  tip- 
over,  and  constant  angle  of  attack  segments  of  the  boost-phase, 
and  assume  that  the  thrust  vectors  of  each  stage  are  applied  at  a 
constant  direction  relative  to  the  inertial  coordinate  system. 
Although  this  is  not  practical  for  an  actual  launching,  this  is  a 
reasonable  approximation  since  the  vertical  and  tip-over  segments 
are  usually  a  small  portion  of  the  boost-phase.  Then  we  assume  that 
the  aerodynamic  drag  is  small  compared  to  the  booster  rocket  thrust, 
and  can  be  neglected.  Finally,  we  neglect  the  variation  of  thrusts 
due  to  altitude,  and  assume  that  the  specific  impulse  and  structural 
factor  are  the  same  for  each  stage. 

The  K- St age  Boost  Socket 

With  these  assumptions  in  mind,  let  us  develop  seme  equations 
for  the  k-stage  rocket  as  shown  in  Fig.  4-2.  Some  useful  relation- 
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T  Uimsr  Vector 
\  Thrust  Any /c. 

Vb  &urnoat  Velocity 

Burnout  Any  It  -  Rect  Coord, 
burnout  A ny/t- Polar  Coord. 
4>w  burnout  Panyt  Arty  fa 


Fi3.4-1 

Coordinate  Cystenr  For  3oosT- 
Phase  Approx/mat/on 


44 


GAW  Mech  61-2 


ships  for  this  rocket  are  listed  as  follows  (Ref  3)*  The 
symbols  used  are  shown  in  Pig.  4-2. 

Structural  factors 

*tll 

initial  mass  of  the  j  stage 
n ,  B  —  -  '  — — 

burnout  mass  of  the  stage 


or 


effective  mass  ratio  of  the 


.th 

0 


sub-rocket  s 


(4-1) 


_  Mj  _  nr\j  +  m}H  j-  ■■■■+■  rnk  -f-  m? _ 

Mjf  -^r  ■+-  mj*'  +  ms+i  +  "  +  "V  ^j~(~n~)  mj 


(4-2) 


where  M  .  is  the  burnout  mass  of  the  j**1  sub-rocket. 

Jf 

stage  mass  ratio: 


(4-j 


overall  mass  ratio: 


r > 


(4— 4a) 


or 


Y.-1+ %  +  *'  +  ■■■+% 


(4-4b) 


Optimization.  It  can  be  shown  that  the  optimum  arrangement  of 
sub— rocket  masses,  which  results  in  a  maximum  final  burnout  velocity 
(vacuum  conditions,  constant  g),  for  a  rocket  of  constant  overall 
mass,  with  the  same  exhaust  velocity  of  each  stage,  is  given  by  the 
equation  (Ref  3) 


0, 


(4-5) 
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If  the  structural  factors  of  each  stage  are  the  sane  then  the  equation 
for  optimization  becomes 


D,  =  D2  =  •  •  •  ~  Dk  (4-6) 

The  overall  mass  ratio  for  optimum  staging,  with  the  stage  structural 
factors  and  exhaust  velocities  equal,  is 


K  -/W 


(4-7) 


The  sub-rocket  mass  ratios  are  then 

%  =  (K - 


(4-8) 


Squations  of  Motion.  In  Chapter  II  it  was  shown  that  the 
equations  of  motion  of  a  one  stage  rocket  in  field-free  space  are 

vk=v0  +  ulnD  (2_7o) 

+  r  (2-76) 

If  we  re-write  these  equations  in  scalar  form  for  our  rectangular 
coordinate  system  and  consider  the  effect  of  gravity  they  become 


xb  =  k0  +  u  cos  X  In  D 

(4-9) 

=  +  u  SLn  Xln  D  -gtb 

(4-10) 

xctk  +  s'  cos  \ 

(4-11) 

y.=  y0tk  +  s  sin  A 

where  s  =0 
o 

=  the  horizontal  component  of  burnout  velocity 
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y^  ■  the  vertical  component  of  burnout  velocity 
xq  =  the  horizontal  component  of  initial  velocity 
yQ  =  the  vertical  component  of  initial  velocity 

X  =  the  angle  between  the  thrust  veotor  and  the  horizon 


and 


(2-76) 


Let  us  now  find  the  components  of  burnout  velocity  and  displacement 
for  a  k-stage  rocket,  starting  from  rest. 


Velocity: 


Stage  1. 

X,  =  u  cos  A  In  D. 

bi 

(4-13) 

ybi  =  u  sin  A  In  D,  -gtfc 

(4-14) 

Stage  2 

kh  =  kh  +  u  cos  A  In  D2 

(4-15) 

Stage  k. 


yb  =  yki  ~h  u  sin  A  In  D2  - gt ^  (4-16) 

X,  =  X,  -h  u  cos  >Zn  0K  (4-17) 


\  =  yv,  -t- u  SLn  * ln  D*  -9% 

(4-18) 

Displacement: 

Stage  1. 

Xb  =  -s'  cos  A 

(4-19) 

,  •  x  sK 

yh)  =  s,  sm  A  -  — ' 

(4-20) 

Stage  2. 

\ = \ + \  \ + K cos  A 

(4-21  i 

-  y»,  +■  y4,  \  +  h  st”  * 

(4-22) 
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Stage  k.  y  =  x.  ~h  X,  t,  -h  s'  COS  A 

kK-<  k*-/  kK  k 


(4-^3) 


>4  ,  +->V  *V,  +-  K  sinb- 


9  < 


’*-/  "  k-l  K 

If  we  now  assume  that  we  have  a  rocket  with  optimum  staging 


(4-24) 


(D's  all  equal),  and  further  assume  that  the  burning  times,  t.  , 

D 

for  all  stages  are  equal ,  we  can  say  that 


y^  =  2s'sin  k  y-^t-,-2^ 


Kbw  from  Bqs  (4-27)  and  (4-28) 

\  +  \  + ' "  +  V 
n ,  +  K -i- -+K. 

But 


\ln  D 

(4-25) 

~2gtk 

(4-26) 

\  In  D 

(4-27) 

-kgt, 

(4-^ 

(4-7  ‘ 

(4-30) 

■t  +"'  +  Xk  )tb 

(4-31) 

(4-32) 

(4-33) 

yh  [l  +  2-h---+(K-l)j 

(4-34  • 

7  _  k/fc-O 
/  2 

(4-35) 
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or 


(4-41) 


•j  ™j 


By  combining  Eqs  (4-40)  and  (4-4l)  we  find  that 


t 


=  ^_  (3tha 

kJ  3  c  ft  Mj 


(4-42) 


But  from  Eq  (4-2)  we  can  show  that 

M,  Dj 

Therefore,  by  combining  Eqs  (4-43)  and  (4-42),  we  find  that 


(4-43) 


t 


i; 


“J  /Dal) 
Sc  &  (  Dj  J 


(4-44) 


We  can  also  show  that 


/  _  Uj  / Dj-l-lnDj ) 
J  9  c  (3j  {  Dj  / 


(4-4,>) 


Prom  Eq  (4-44)  it  is  seen  that  if  the  exhaust  velocities,  sub-rocket 

effective  mass  ratios,  and  thrust  to  initial  weight  ratios  are  all  the 

same,  then  the  burning  times  are  also  all  the  same.  This  is 

exactly  the  assumption  we  made  to  get  Eqs  (4-27 »  28,  36  and  37).  By 

using  Eqs  (4-44)  and  (4-45)  it  can  be  shown  that  the  final  velocity 

th 

and  displacements  at  burnout  of  the  k  stage  are 

X.  =  kuln  D  cos  X  (4-46) 

6k 

yk^  =  ku  In  D  sin  X  -  (4-47) 

ku*cosX_  r ZCD'l-ln  D)  +  (k-0(P-Qln  D  1 
£S*/S  j  D  J  (4-48) 
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kuz  sin  A  f  2(D-t-ln  D)  +  (k-0(D-!)ln  0 

29c/£?  7  D  J 


/  r ku(D-l)  T 
29c[  P  D  J  (4-49) 

Thus,  with  the  burning  time  eliminated  from  the  equations  hy  use 
of  the  thrust  to  weight  ratio,  we  oan  say  that 

Burnout  Conditions  «  f(k,  u,  X,  j5,  D)  (4-50) 

We  must  now  combine  these  equations  with  the  correoted  trajectory 
equations  of  Chapter  III,  so  that  we  may  find  a  method  for  evaluating 
a  function  of  the  form 

I)  «  f(k,  u,  X,  p,  Burnout  Conditions)  ( 4—51 ) 


Technique  for  Satisfying  Burnout  Velocity  Requirement 

First  let  us  solve  the  problem  hy  satisfying  only  the  conditions 
of  burnout  velocity  and  direction,  and  later  consider  what  must  be 
done  to  also  satisfy  the  condition  of  burnout  position.  Assume  that 
we  are  given  a  required  burnout  velocity,  v^,  and  burnout  angle, 
as  shown  in  Fig.  4-1 »  and  that  we  want  to  find  the  required  sub-rocket 
mass  ratios,  D.  We  assume  also  that  k,  and  £3  are  given.  We  can  say 

that 


tan  77 
*  y 


(4-52; 


where  and  y^  are  given  by  Eqns  (4-46)  and  (4-47)  with  the 
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k  subscript s  omitted  for  simplicity.  Then 

tan  Tjt  =  tan  A  - 

If  we  multiply  both  sides  of  Eq  (4-53)  by  the  quantity,  cos  A  oos  r) 
and  rearrange ,  the  resulting  equation  is 

\  \  (D-l)  cos  il* 

stn  A  cos  7fk  -  cos  A  sin  7^  -  —^l^D —  (4-54) 


Then 

A  \  (D-l)  cos  7?a 

sin(A-Vk 

(4-5 4a) 

Solving  Eq  (4-54a)  for  A, 

V  •  -ir(D-l)  cos  TUI  . 

sm  L  fiD  VTd~J  +  X 

(4-^>) 

or 

\  -  ??,  +■  sin''  f 

(4-5:) 

where 

K>  (D-0  COS  Tlk 

S  fiD  In  D 

(4-57) 

From  Fig.  4-1  the  horizontal  component  of  burnout  velocity  is 

h  =  Vh  cos  Vk  (4-58) 

By  combining  Eqs  (4-46,  55  >  57 ,  and  58)  we  can  show  that 

ku  In  D  cos  (yk-h  sin'1  f)  -  vk  cos  n lh  =■  O  (4-59) 


or  -F(D)=  O 

Now  we  have  an  equation  with  one  unknown,  D,  which  would  be 


(4-60) 
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difficult,  if  not  impossible,  to  solve  explicitly  for  D.  A  numerical 
solution  is  found  by  using  Newton's  iteration  method  (Ref  14* 495 ) • 

With  this  method,  we  choose  a  trial  value  for  E,  say  For  this 

value  we  evaluate  f(B^)  and  f'^)  where 


The  equation  for  f'(B)  can  be  shown  to  be 

rrD)*ku^ 


( D-l)  cos  A  f  (V7  ~  In  D)  sin  A 

D(D- 1)  — 


(4-6la) 


Then  our  new  trial  value  for  H  is  given  by  the  formula 


£ (Dj) 
■f'(Oi) 


(4-62) 


This  process  converges  quite  rapidly  in  this  particular  case,  and 
is  repeated  until  the  required  degree  of  accuracy  is  attained.  Thus 
we  can  find  D,  the  sub-rocket  effective  mass  ratio  required  to 
achieve  a  given  burnout  velocity  and  direction. 


Technique  for  Satisfying  Burnout  Position  Requirement 

In  order  to  simultaneously  satisfy  the  oondition  of  burnout 

altitude,  y^,  or  radius,  r^,  we  must  resort  to  a  second  trial  and 

error  process.  To  do  this,  first  choose  a  trial  value  for  burnout 

range  angle,  9^,  (See  Fig.  4-3a).  With  the  chosen  4>  ^  and  the 

corrected  interceptor  trajectory  parameters,  we  use  the  general 

trajectory  equations  of  Chapter  II  to  determine  rv  ,  v,  ,  and  v,  . 

°1  °1  °1 

As  a  first  approximation  let  us  assume  that  the  launch  site  in  the 
rectangular  coordinate  system  coincides  with  the  theoretical  launch 
site  in  the  polar  coordinate  system.  Then  we  can  say  that 
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^nioa 


(4-63) 


where. these  angles  are  shown  in  Fig,  4-3a.  Therefore  the  first 
trial  hurnout  angle  in  rectangular  coordinates  is 

^  =  %-€,  (4"64) 

Then  by  using  Eq  (4-59)  we  can  find  a  value  for  With  known, 

Eqs  (4-48)  and  (4-49)  are  used  to  find  x,  and  y,  . 

°1  bl 

We  can  now  compute  the  actual  value  of  r^,  which  is  designated  r^  . 
using  the  relationship  (See  Fig.  4-3a) 


=  V^+yJ 


*  +  4* 


(4-65) 


Now  we  can  remove  the  restriction  that  the  rectangular  and  polar 
launch  sites  are  coincident  by  finding  that 


=  tan 


-I 


rn  +■ 


A 


(4-66) 


and  change  from  the  coordinate  system  of  Fig.  4-3a  to  the  adjusted 

coordinate  system  of  Fig.  4-3b.  With  the  value  of  r.  from  Eq  (4-65) 

°2 

and  the  interceptor  trajectory  equations  we  can  find  new  values  of 
4> ,  ,  v,  .  and  v.  .  Now  substituting  €.  and  v,  into  Eq  (4-64)  we 

b2  V  *2  2  *2  ’ 

find  ti^  .  Then  again  using  Eqs  (4-59?  48,  and  49)  we  find  values 

for  Dg,  and  y^  .  The  iteration  process  is  continued  until 


(4-67) 


Note  that  we  use  the  previous  value  of  6  to  compute  a  new  value 
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for  rj  .  This  introduces  a  small  error  for  the  first  few  trials 
b 

since  in  all  likelihood  the  relative  positions  of  the  rectangular 

and,  polar  launch  sites  will  be  different  from  one  trial  to  the 

next.  However,  since  this  process  is  convergent,  as  the  number 

of  trials  increases,  €  will  approach  a  constant  value  and  the 

error  will  approach  zero.  An  exaggerated  illustration  of  how  this 

iteration  process  might  converge  is  shown  in  Pig.  4-4*  Thus,  by- 

using  this  double  iteration  process,  we  find  the  effective  sub-rocket 

mass  ratio,  D,  for  the  boost-phase  rockets,  and  values  for  4>  and  r,  . 

b  b 

We  are  now  prepared  to  achieve  our  objective,  and  determine 
the  weights  and  thrusts  required,  and  calculate  the  excess  time. 

First  let  us  find  the  weights  and  thrusts. 


Determination  of  Weight  and  Thrust  Requirements 

The  effective  sub-rocket  mass  ratio  for  the  velocity-match 
rocket  is  found  by  using  Eq  (3-38).  Let  us  also  specify  a  thrust 
to  initial  weight  ratio  for  this  rocket,  and  call  this  quantity 
Pu,  where  the  u  subscript  is  used  to  avoid  confusion  with  the 
boost-rocket  parameters.  Assuming  a  one-stage  velocity-match 
rocket,  the  overall  mass  ratio  of  this  rocket  is  found  from  Eq  (4-7) 
and  Pig.  4-2  to  be 


m'p  _  (n-Q  Du 

(n-Du) 


(4-68) 


where  ^  is  the  non-dimensional  mass  (or  weight)  ratio  of  the 
velocity-match  rocket  (lb  of  int.  rocket  per  lb  of  int.  vehicle). 
The  thrust  required  for  the  velocity-match  rocket  is  found  from 
Eq  (4-38)  to  be 


=  ft 


(4-69) 
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or  in  non-dimensional  form 

Thu  -  pu  %  (4-70) 

where  Ty  -  the  thrust  of  the  velocity-match  rocket 

-  the  weight  (or  lb  mass)  of  the  velocity-match  rocket 
Th^  «  the  lb  of  thrust  required  per  lb  of  interceptor  vehicle. 
The  mass  ratios  of  the  boost-stages  are  found  by  using  Eqs  (4-7) 
and  (4-8).  These  mass  ratios  are 


(4-7) 

(4-71) 

-  (r°  wf1  vir. J 

(4-72) 

<.4-73) 

The  thrusts  required  for  each  stage  are,  from  Eq  (4-38)  and  Fig.  4-2 


T,  =  (3  M,  =  p(m,  +  m2+'-  +  mk+m'f) 

(4-74) 

@ (m2  -h  m3  -f —  +  mk  j- 

(4-75) 

Tk=  ?Mk=  p(mK  +  m,) 

(4-76) 

By  dividing  through  by  m\  and  using  Eq  (4-68)  we  can  write 
these  equations  in  non-dimensional  form. 
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Th,  -  &  r„(r,  +  %  +  ■■■  +  rK  + 1) 

(4-77) 

Th,  =  /sru(%-h%  +  -+rK  +  i) 

(4-78) 

77 =  + 

(4-79) 

By  using  Eqs  (4— 4b )  and  (4-8)  we  can  further  reduce  these  equations 
and  show  that 


Th ,  = 

(4-80) 

-P,  _ 

(4-81) 

T^=/f  k  t  j 

(4-82) 

Therefore,  the  mass  (or  weight)  and  thrust  ratios  are  determined, 
and  numerical  values  of  weights  and  thrusts  can  he  found  if  m^, 
the  weight  of  the  interceptor  vehicle,  and  k,  the  number  of 
boost-phase  stages  are  known. 

Excess  Time 

To  find  the  excess  time  it  is  necessary  to  determine  the  flight 
times  to  the  intercept  point  for  both  the  warhead  and  the  interceptor 
The  difference  between  these  times  is  defined  as  the  excess  time. 

Let  us  assume  that  the  burning  time  of  the  enemy  booster 
rockets  can  be  neglected  and  that  the  total  flight  time  of  the 
warhead  from  its  launch  site  to  the  intercept  point  is  the  same  as 
its  free-flight  time.  This  time  can  be  found  by  using  the  warhead 
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trajectory  parameters  and  evaluating  Eqs  (2-46)  and  (2-47)  at  the 

* 

angles  eh  »  and  <h  as  shown  in  Fig.  4-5*  Let  this  time  he  T  . 

The  time  of  the  interceptor  from  its  launch  site,  LS^,  to  the 

boost-phase  burnout  point,  BO,  can  be  found  by  multiplying  Eq  (4-44) 

* 

by  K,  the  number  of  booster  stages.  Let  this  time  be  T^  .  From 

Fig.  3-5  it  can  be  seen  that  the  flight-time  of  the  interceptor 

from  point  A'  to  the  intercept  point,  IP,  is  the  same  as  the  free- 

flight  time  (without  velocity-match  rocket  ignited)  from  point  A* 

to  the  aiming  point,  AP.  Thus  the  free-flight  time  from  point 

BO  to  IP  (Fig.  4-5)  can  be  found  by  using  the  intercept  trajectory 

parameters  and  evaluating  Eqs  (2-46)  and  (2-47)  at  the  angles  ep  , 

*  * 

and  4>  .  Let  this  time  be  T  .  Then  the  excess  time,  T  , 

r  inxa  c  x 

is  given  by  the  equation 


*  ,  *  *  x 
T  -  (T,  +  T  ) 

w  b  c 


(4-83) 


Therefore,  by  assuming  that  the  enemy  missile  is  detected  at  the  instant 

* 

it  is  launched,  the  interceptor  must  be  launched  T^  seconds  later,  if  the 
warhead  and  interceptor  vehicle  are  to  arrive  at  the  intercept  point 
simultaneously. 


Maximum  Acceleration 

An  important  consideration,  especially  if  a  manned  vehicle  is  used 
in  the  interceptor  system,  is  the  maximum  acceleration  imposed  upon  the 
system.  There  is  a  peak  acceleration  at  the  instant  of  burnout  of 
each  sub-rocket,  and  this  is  found  by  using  Newton's  second  law  to  be 
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where 


v  . 
J 


T. 

J 


the  maximum  acceleration  in  sea  level  "g's" 

tb. 

the  thrust  of  the  j  stage 

"til 

the  burnout  weight  of  the  j  stage. 


But.  from  Eq  (4-38) 


Tj  -  /Sy  Mj 


(4-85) 


By  combining  Eqs  (4-84)  and  (4-85)  it  is  found  that 

4-  Mj 


But 


My 


M 


'--Dy 


By  combining-  Eq  (4-2)  ar.d  (4-86)  we  get 

ft  =  ft  ft 


(4-86) 


(4-2) 


(4-87) 


Since  we  have  specified  that-  the  ?’s  and  P’s  of  each  sub-rocket 
are  equal,  Eq  (4—87)  shows  us  that  the  peel:  accelerations  during 
the  boost-phase  are  the  same.  Similarly,  the  maximum  acceleration 
during  the  velocity-match  phase  would  then  be 


\>U  -  Du 

During  an  actual  launching  there  would  be  other  components  of 
acceleration  (eg.  angular  and  centripetal  accelerations  during  the 
tip-over  maneuver)  but  these  are  considered  to  be  of  secondary 
importance  for  this  investigation. 


Summary 

In  this  chapter  a  method  for  determining  the  booster-rocket 
requirements  was  developed.  This  was  accomplished  by  making  several 
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simplifying  assumptions,  developing  the  equations  of  motion  of  a 
k-stage  rocket  and  introducing  the  concept  of  thrust  to  weight 
ratio.  The  equations  derived  in  this  manner  were  then  combined 
with  the  interceptor  trajectory  equations  of  Chapter  III  to 
determine  thrust  and  weight  ratios,  available  warning  time,  and 
maximum  accelerations  of  the  boost-phase  and  velocity-match 
rockets. 

Because  of  the  many  equations,  and  two  iteration  processes 
involved,  numerical  results  are  very  difficult  to  obtain  without 
the  aid  of  an  electronic  computer.  For  this  reason,  the  equations 
are  programmed  on  the  IBM  1620  digital  computer  to  permit  the 
investigation  of  many  different  intercept  conditions.  The 
computer  programs  used  in  the  analysis  are  described  in  Appendix  A. 
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V.  Procedure  for  Obtaining  Numerical  Results 

Primary  Considerations 

The  overall  vehicle  weight,  the  thrust  of  the  first  boost  stage, 
the  excess  time,  and  the  maximum  "g"  forces  on  the  interceptor  vehicle* 
are  the  principal  factors  to  he  determined. 

The  excess  time  is  certainly  the  most  important  factor  in  the 
analysis,  for  without  adequate  excess  time,  the  intercept  technique 
is  impractical  or  even  impossible.  Therefore  we  wish  to  determine 
what  intercept  conditions  allow  the  most  excess  time,  and  the  effects 
of  the  boost  parameters  on  excess  time. 

The  overall  vehicle  weight  is  also  quite  important,  since  as 
weight  increases,  the  thrust  xoquiredto  launch  the  vehicle  increases. 
Since  there  are  limits  on  the  amount  of  thrust  available  to  perform 
the  intercept,  thrust  rather  than  weight  is  the  factor  which  we  atteiiq 
to  minimize  in  the  analysis. 

Finally,  the  accelerations  which  exist  during  the  boost  phases 
are  also  an  important  factor  if  the  system  is  manned.  Man  is 
capable  of  withstanding  approximately  ten  transverse  "g’s"  or  322 
feet  per  second  per  second  applied  with  the  pilot  perpendicular  to 
tne  acceleration  vector  (Ref.  4:14).  Therefore  we  shall  indicate 
the  intercept  conditions  which  result  in  accelerations  in  excess  of 
ten  *'g’ s'*. 

R*3ic  Analysis 

The  basic  analysis  consists  of  fixing  the  warhead  trajectory, 
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intercept  point,  and  the  boost  parameters  and  then  varying  the 
position  of  the  interceptor  launch  site.  This  procedure  is 
repeated  for  different  intercept  points  on  the  warhead  trajectory. 
Finally  the  complete  procedure  is  repeated  for  other  warhead 
trajectories.  Once  the  basic  analysis  is  complete,  the  effects  of 
varying  I  with  n  fixed,  varying  n  with  I  fixed,  and  varying  I 

S  8  S 

and  n  in  combination  are  determined. 
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It  was  shown  in  Chapter  II,  that  for  any  given  range  there 
“is  an  almost  unlimited  number  of  trajectories,  depending  upon  the 
amount  of  total  energy  available.  Thus  we  must  stipulate  more  than, 
just  the  range  to  define  the  exact  trajectory.  In  order  to  more 
completely  consider  the  possible  intercept  situations,  three  different 
trajectories  are  chosen  for  each  range.  One  is  the  minimum-energy 
trajectory,  and  the  others  are  a  high  angle  (high  altitude)  and  a  low 
angle  (low  altitude)  trajectory.  These  trajectories  can  be  completely 
defined  if  we  now  specify  the  theoretical  launch  angle,  y0>  of  the 
warhead.  For  the  minimum  energy  trajectories,  we  can  use  Eq.  (2-57) 
to  determine  yo*  For  the  6000  mile  trajectory 

tan  Yq  -  C08<fr°  *  0.64278759  =  0.36192092 

m  1  +  sin  4>o  1.76604431 

and  y  *  0.34906588  radian. 

°ME 

For  the  4000  mile  trajectory, 

tan  y  »  °-8354S776  ,  0.53919520 

°ME  1.54950890 

and  Ya  *  0.49450997  radian. 

TiE 

To  establish  the  high  and  low  angle  warhead  trajectories  for 
each  range,  we  can  fix  the  total  energy  available  for  the  trajectory. 
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This  than  determines  the  energy  ratio  of  the  trajectory,  which  in  tura 
determines  the  theoretical  launch  angle.  Since  the  choice  of  total 
energy  is  an  arbitrary  one  as  long  as  reasonable  values  are  used,  it 
is  more  convenient  to  specify  the  energy  ratio.  An  energy  ratio  of 
0,9  is  chosen  for  the  6000  mile  trajectory,  and  0.8  for  the  4000  mile 
trajectory.  Then  for  the  6000  mile  trajectory 


tan  T  Q  = 


1 

2 


ER  ^ 

if  m  y  -  4(  1  -  er  ) 

tan  0 

Vi  tan 

0.9 

t] 

1 

h 

3 

• 

CXI 

ON 

• 

0 

0.9  ) 

1.1917534 

((1.1917534/ 

(2-58) 


0.58389717  and  0.17126302 


or  yo  =  0.52849701  radian  (high  angle) 

and  "if  0  =  O.I6961984  radian  (low  angle) 


For  the  4000  mile  trajectory 


1 

2 


0.8 


0.8 


0.65771029  yl  0.65771029 ) 


-  4(  1  -  0.8  ) 


«  I.02032565  and  O.I96OI585 

or  if  =  0.79544624  radian  (high  angle) 

o 

and  -g-  *  O.I9356362  radian  (low  angle) 
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These  values  of  4>  and  y  which  are  listed  in  Table  I, 
are  necessary  inputs  to  the  computer  solution  of  the  boost-phase 
approximations  (Appendix  A),  since  the  computer  is  programmed  to 
determine  all  other  parameters  of  the  warhead  trajectory  from 
these  values. 


Table  I 

Warhead  Trajectory  Semi-Range  and  Launch  Angles 
for  6000  and  4000  Nautical  Mile  Ranges 


Warhead 
Range 
(n.m. ) 

Trajectory 

Angle 

Semi-Range 

Angle 

(radians) 

Launch 

Angle  ' 

(radians) 

Low 

0.58177638 

0.19356362 

4000 

Minimum  Energy 

It 

0.49450997 

High 

II 

0.79544624 

Low 

0.87266456 

O.I6961984 

6000 

Minimum  Energy 

M 

0.34906588 

High 

tl 

0.52849701 

Selection  of  Intercept  Points  and  Launch  Sites 

To  facilitate  the  discussion  of  the  position  of  intercept 
points  and  launch  sites,  we  will  use  the  angular  position  of  the 
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points  relative  to  the  warhead  apogee  in  an  earth-centered,  polar 
coordinate  system.  This  is  illustrated  in  Pig.  5-1  • 


Polar  Axis 


Fi  g.  S-  I 

Position  of  Launch  Site  an<i 
Tnt ercept  Point  Relative  to 
Warhead  Apogee 


An  intercept  point  is  defined  by  <£  and  is  understood 

to  be  on  the  warhead  trajectory.  The  position  of  the  interceptor 

launch  site,  </>  T  <,  ,  is  determined  in  the  same  manner,  but  its  position 
ra 

vs  on  the  earth's  surface.  Note  that  t>  is  positive  clockwise 
and  the  warhead  apogee  is  at  -  0.  With  these  definitions 

established  we  can  proceed  with  the  selection  of  intercept  points  and  launch 
sites. 


The  intercept  points  to  be  investigated  are  arbitrarily  chosen 
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so  that  a  represent at.i vs  sample  of  points  on  the  warhead  trajectory 

is  cheeked.  These  points  are  approximately  equally  spaced  along 

each  trajectory  between  warhead  apogee  and  atmospheric  re-entry  point. 

The  interceptor  launch  sites  are  determined  indirectly  by 

using  the  theoretical  intercept  trajectory  semi-range  angles  as 

inputs  to  the  computer  program.  Then  the  approximate  position 

of  the  interceptor  launch  site  relative  to  the  warhead  apogee  is 

S  -  <b  .  The  exact  position,  which  is  4>  -  <d> 

Y  wr  r  so  ~wx  r* 

is  slightly  further  up-range  than  the  approximate  position  (Fig.  3-6) > 
hut  this  is  not  critical  since  we  are  merely  trying  to  select  a 
representative  group  of  launch  site  positions.  The  up-range 
positions  of  interceptor  launch  sites  are  limited  by  practical 
considerations;  e.g.  no  launch  sites  in  enemy  territory.  The 
down-range  position  of  the  interceptor  launch  site  is  also  limited. 

As  the  angular  position  of  the  launch  site  approaches  the  intercept 
point,  the  interceptor  trajectory  becomes  almost  vertical  and  the 
trajectory-match  technique  is  not  practical. 

With  these  limitations  in  mind,  the  interceptor  launch  sites 
are  arbitrarily  chosen  to  determine  enough  points  to  plot  curves 
which  show  the  effect  of  launch  site  position  on  the  boost  requirements. 

Selection  of  Rocket  Parameters 

The  five  remaining  variables  to  be  investigated  are: 

(1)  k,  the  number  of  stages  for  the  boost  rocket 

(2)  I  ,  the  specific  impulse  of  the  fuel 

s 

(3)  n,  the  structural  factor 
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(4)  P ,  'the  thrust  to  weight  ratio  for  the  boost-phase 

(5)  P  ,  the  thrust  to  weight  ratio  for  the  velocity-match 

rocket . 

In  order  to  continue  the  investigation  some  particular  values  of 

these  variables  must  be  selected. 

Number  of  Boost-Phase  Stages.  For  this  investigation  three 

boost  stages  were  selected  as  a  compromise  between  the  decreased 

reliability  of  a  large  number  of  stages  and  the  large  thrust  and 

weight  requirements  for  a  small  number  of  stages.  This  reduces  the 

number  of  variable  boost  parameters  to  four. 

Specific  Impulse  and  Structural  Factor.  Specific  impulse  is 

determined  by  the  choice  of  propellant.  Values  of  I  vary  from 

s 

160  seconds  for  some  old-type  solid  propellants,  to  364  seconds 
for  liquid  hydrogen-liquid  oxygen  (Ref  10:  112,  113,  312,  313)* 
Instead  of  selecting  a  particular  fuel,  a  reasonable  value  of  I  , 

9 

2S0  seconds,  is  selected  for  the  basic  analysis.  Later  in  the 
investigation  the  effect  of  different  values  of  Ig  for  a  few 
typical  intercept  situations  is  determined. 

A  similar  approach  is  used  to  determine  the  value  of  the 
structural  factor  since  it  depends  upon  the  density  of  the  fuel 
and  the  detailed  structural  design  of  a  rocket.  A  value  of  n  equal 
to  xen  is  selected  for  the  basic  analysis. 

Thrust  to  Initial  'weight  Ratios.  The  choice  of  numerical 
values  for  p  and  appears  to  be  completely  arbitrary.  Then  to 
determine  the  values  to  be  used  in  the  analysis,  it  is  necessary  to 
fix  the  value  for  one  at  an  arbitrary  value,  and  determine  the  effect 
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of  the  other  on  overall  weight  and  thrust  requirements.  Once  thes- 
effects  are  known,  the  most  promising  values  are  used  in  the  "basic 
analysis. 


Summai 


This  chapter  outlines  the  general  approach  used  to  obtain  the 


desired  data  from  the  digital  computer  solution  of  the  intercept 


analysis.  These  results  are  presented  in  Chapter  VI. 
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VI.  Results  and  Conclusions 


A  numerical  integration  was  performed  using  computer  program 
CM-12  (Appendix  A),  for  one  particular  intercept  trajeotory.  In 
this  integration,  the  effects  of  drag  and  the  vertical,  constant 
angle  of  attack,  zero  lift,  and  constant  attitude  flight  segments 
of  a  rocket  launching  are  considered.  This  was  accomplished  to 
verify  the  validity  of  the  assumptions  made  in  the  "boost-phase 
analysis.  For  the  case  considered,  the  actual  burnout  velocity 
came  within  three  percent  of  that  required.  The  results  of  this 
integration  are  summarized  in  Appendix  B. 

Since  tne  results  are  only  good  approximations  to  the  actual 
’-alues  because  of  the  many  simplifying  assumptions  made,  the  coaputt: 
data  is  presented  in  graphical  rather  than  tabular  form.  This  also 
gives  a  better  overall  picture  of  this  intercept  technique,  and 
allows  important  trends  to  be  interpreted  more  easily.  Angular 
inputs  and  print-outs  of  the  computer  are  in  radians  and  therefore 
all  angles  are  plotted  on  the  graphs  in  radians. 

In  order  to  avoid  complicating  some  of  the  graphs,  a  code 
is  used  to  identify  the  intercept  conditions.  For  example,  6H$45s30 
la  interpreted  as  follows: 

6H  6000  nautical  mile,  high  angle  warhead  trajectory 

45  intercept  point  on  warhead  trajectory  is  0.45  radians 

down  range  from  the  apogee 

30  3emi-range  angle  of  the  theoretical  intercept  trajectory 
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is  0.30  radians 

The  theoretical  semi-range  angle  is  selected  for  this  code  rather 
than  the  corrected  value,  because  this  is  the  angle  used  a3  the 
computer  input  to  obtain  the  numerical  results  which  are  plotted. 

Results 

Effect  of  on  Th^.  Figures  6-1,  2,  3  were  obtained  for 
two  intercept  conditions:  6H :45s 30  and  6H: 45*60.  Thrust  to 
interceptor  vehicle  weight  ratio,  p  ,  was  varied  while  the  other 
rocket  parameters  were  held  constant.  These  values  are 

I  ■  250  seconds 

s 

n  =  10 

3  »  1.5 

k  =  3 

For  the  case,  6H: 45*60,  has  a  negligible  effect  on  Th^  (Fig.  6-1 ). 

thus  little  advantage  is  gained  by  using  a  more  powerful  engine  for 

the  velocity-match  rocket.  The  top  curve  shows  that  a  reduction 

in  thrust,  Th^,  of  approximately  3*5^  can  be  obtained  by  increasing 

from  one  to  three;  however,  since  the  curve  appears  to  be 

leveling  off,  little  will  be  gained  by  using  a  more  powerful  engine. 

A  comparison  of  the  curves  in  Fig.  6-1  shows  that  p^  has  more  effect 

for  the  intercept  trajectory  with  the  smaller  theoretical  semi-rangf 

angle  (.3  vs  .6).  This  is  due  to  the  larger  difference,  Av  ,  between 

warhead  and  interceptor  velocities.  A  larger  decreases  the 

distance  s'  (See  Eq  4-45  and  Fig.  3-6).  This  also  decreases  rT  an 

xa. 

therefore  we  can  say  that 
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lim  (corrected  int.  trajectory)  =--  (theoretical  int.  trajectory) 
(3U  — *  ao 

In  other  words  we  can  say  that  as  r  approaches  rT ,  the  total  energy, 

la  J- 

E,  of  the  intercept  trajectory  approaches  a  constant  value  (E  of  the 
theoretical  intercept  trajectory),  thus  explaining  why  the  curves 
in  Fig.  6-1  level  off. 

Effect  of  p  and  |3  on  v.  The  computer  data  indicated  that  for 
the  values  of  |3  used,  the  maximum  acceleration  during  the  boost-phase 
ranged  from  approximately  3*0  to  4*5  "g’s",  and  thus  this  acceleration 
is  not  considered  to  be  a  problem.  Therefore  no  graphs  are  plotted 
for  this  data. 

Figure  6-2  shows  that  the  maximum  acceleration  of  the  velocity- 
match  rocket  is  a  linear  function  of  p  .  This  is  also  apparent 
from  E^  (4-38), 

^ u~  @u 

where  is  a  constant  fox'  a  particular  intercept  point.  Since  Du 
is  a  function  of  Av,  and  Av  increases  with  decreasing  semi-range 
angles,  the  maximum  acceleration  can  be  a  limiting  factor  on  a 
choice  of  for  intercept  trajectories  with  small  semi-range  angles, 
especially  for  a  manned  interceptor  vehicle.  For  the  two  cases 
shown  in  Fig.  6-2,  it  can  be  seen  that  for  a  given  the 
acceleration  is  more  than  doubled,  when  the  serai-range  angle  is 
halved. 

Effect  of  on  T  .  The  same  reasoning  that  was  used  to 

* 

analyze  the  effect  of  on  Th^  is  used  for  the  effect  of  on  T^ 
(Fig.  6-3).  Here  again,  the  effect  of  j3^  is  more  pronounced  for 
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the  smaller  semi -range  angle,  and  both  curves  will  eventually 
approach  a  maximum  value  of  excess  time.  An  inorease  in  from 
one  to  three  results  in  an  increase  in  T  of  only  about  3 %  and  2.5  $ 
for  the  upper  and  lower  curves  respectively.  Although  this  small 
percent  change  might  be  important  in  some  cases,  we  can  safely  say 
that  in  general  Pu  is  not  an  important  factor  as  far  as  time  is 
concerned. 

Effect  of  p  on  Th^.  Figures  6-4  and  6-5  are  plotted  from 
computer  data  obtained  by  varying  p  and  holding  the  other  variables 
fixed  as  follows: 

I  -  250  seconds 

s 

n  -  10 

P  =  3.0 

u 

k  «  3 

Figure  6-4  shows  that  for  each  set  of  intercept  conditions  there  is 
an  optimum  value  of  p  which  results  in  a  minimum  value  of  Th^. 

These  minima  are  1.8  and  1.6  for  the  upper  and  lower  curves 
respectively.  However,  since  the  curves  are  quite  flat  in  these 
regions,  a  value  of  1.5  was  arbitrarily  chosen  as  a  representative 
value  of  the  optimum  p.  The  existence  of  these  optimum  values  can  be 
explained  by  considering  Eq  (4-44),  which  shows  that  the  burning 
times  of  the  boost-phase  stages  are  inversely  proportional  to  p, 
and  the  equation 


Th ,  =  p  %  % 


(4-80) 
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As  p  is  first  increased,  the  burning  time  shortens,  and  the 

velocity  loss  due  to  gravity  (  gt^  )  decreases,  thus  less  fuel 

is  required  to  reach  a  given  velocity  and  decreases.  Hov/ever, 

a  point  is  reached  where  p  increases  faster  than  decreases.  At 

this  point  the  thrust  required,  Thn ,  is  a  minimum. 

* 

Effect  of  p  on  T  .  Figures  6-5  shows  that  an  increase  in  p 
results  in  a  larger  excess  time,  and  that  this  time  approaches  some 
limit  asymptotically.  The  same  explanation  that  was  used  in 
discussing  the  effect  of  P^  applies  here  also.  We  can  say  that 

lim  (interceptor  flight  time)  =  (free  flight  trajectory  time) 
ft  — >  oo 

which  means  that  for  an  infinite  acceleration  we  can  attain  the 

theoretical  launch  velocity  at  the  earth's  surface  and  "coast" 

along  the  trajectory  from  the  ground  up.  An  increase  in  p  from  1.2 

to  1.5  results  in  increases  in  of  approximately  10$  and  28$  for 

the  upper  and  lower  curves  respectively.  This  means  that  a 

considerable  advantage  is  gained  by  using  a  high  value  of  p ,  for  the 

intercept  trajectories  with  larger  semi-range  angles,  ,  where 

* 

T^  might  be  critical. 

Effect  of  n  on  Th^.  The  data  used  to  plot  the  graphs  of  Th^  vs 
n  (Fig.  6-6)  was  obtained  by  holding  the  other  variables  fixed  as 
follows: 

I  =  250  seconds 

s 

Pu  =  3.0 

p  =  1.5 
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1.0  1.2  1.4  1.6  1.8  2.0 


Thrust  to  Initial  Weight  Ratio  of 
Boost-Phase  Sub-Rockets  (3 

Fig.  6-5 

_  * 

The  Effect  of  (3  on  T 
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k  -  3 

The  curves  show  that  advanced  structural  design  techniques  (higher  n) , 
result  in  a  great  reduction  in  Th^,  hut  that  a  point  of  diminishing 
returns  will  he  reached,  and  the  curves  approach  a  limiting'  value. 

This  limiting  value  is  unattainable  since  as  n  approaches  infinity, 
the  stage  weight  is  equal  to  the  weight  of  fuel.  This  assumes  that 
no  structure  is  required  to  carry  this  fuel.  The  structural  factor 
has  no  effect  whatsoever  on  excess  time. 

Effect  of  I  on  Th  .  The  curves  in  Figs.  6-7  and  6-8  were 

8  i- 

obtained  from  computer  data  with  the  other  variables  fixed  as  follows: 

eu  «  3.o 

P  -  1.5 

n  *  10 

K  =  3 

Figure  6-7  clearly  shows  why  great  emphasis  is  placed  upon  discovering 
new  and  exotic  rocket  fuels.  For  example,  an  increase  in  I  from  200 
to  225  seconds  results  in  decreased  thrust  requirements  of  approximately 

5C$,  and  an  increase  in  I  from  225  "to  250  seconds  results  in  a  further 

8 

decrease  of  approximately  40 ‘fo. 

Effect  of  I  on  T  .  Figure  6-8  shows  that  T  is  practically 

8  X  X 

a  linear  function  of  I  .  This  can  also  be  seen  from  the  equation 

s 

f  _ — yi— 

b  3iP\  D  J  (4-44) 

or  -a  _  f  D-  1  1 

e  (  d  J 
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Pi«.  6-8 

The  Sffeet  of  I#  on 
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Thus  as  Ig  increases,  burning  time  increases  which  results  in  a 
greater  velocity  loss  due  to  gravity,  and  therefore  causes  a 
decrease  in  T  .  For  the  intercept  condition  6H:45J60,  an  increase 
in  I  from  200  to  225  seconds  decreases  T  by  approximately  2.5$> 

8  X 

but  this  small  change  is  insignificant  when  compared  to  the  50$ 
saving  in  Th^. 

Combined  Effect  of  Changes  in  I  and  n.  For  comparison 

purposes,  the  effect  of  combinations  of  optimistic  and  pessimistic  values 

of  I  and  n  upon  Th,  for  the  6H: 45*60  intercept  are  determined, 
s  1 

These  results  are  shown  in  Table  II. 


Table  II 

Combined  Effects  of 

Specific  Impulse  and  Structural  Factor 
on  Overall  Thrust  Requirements"' 


Specific 

Structural 

Thl 

Impulse  (sec.) 

Fact  or 

8 

796 

200 

20 

251 

8 

45 

350 

20 

30 

*  Note:  For  6H:45s60  intercept,  (3  =  1.5.  =  3.0 
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Effect  of  <fi  LS  and  <f>  on  Th^  and  T  .  Figures  6-9,  10, 
m 


and  11  show  the  effect  of  launch  site  position,  TC(  ,  on  Th,  for 

Lb  1 

m 

various  intercept  points  on  the  6000  nautical  mile  warhead  trajectories. 

*  , 

Figures  6-12,  13  and  14  are  the  curves  of  T  vs  9TO  which  correspond 

x  hbm 

to  the  Th^  curves.  On  the  Th^  graphs,  the  dotted  lines  ( . 


) 


indicate  t-iie  approximate  position  of  lines  of  constant  T  (isochrones )  = 
The  points  marked  with  a  star  (  *  )  are  those  points  where  the 
maximum  acceleration  of  the  velocity-match  rocket  reached  10  "g's". 

For  all  points  above  the  star  on  each  curve,  is  greater  than  10 
"g's". 

Many  of  the  Th^  curves  appear  to  approach  horizontal  and  vertical 

assymptotes.  As  the  launch  site  is  moved  up-range  from  the  warhead 

apogee,  mQ&  approaches  41  the  interceptor  and  warhead 

trajectories  will  eventually  coincide  and  determine  the  horizontal 

limit  for  Th, .  The  vertical  limit  occurs  when  4>  is  equal  to  zero. 

1  '  mo 

However,  this  is  simply  a  vertical  trajectory.  The  interceptor 

velocity,  v  ,  at  the  intercept  point  therefore  approaches  zero  as 

<^>  approaches  zero,  which  results  in  a  high  Av,  and  consequently  a 

much  larger  velocity-match  rocket  is  required.  Thus  a  much  heavier 

boost-phase  rocket  with  higher  thrust  is  required  for  these  cases. 

Illustrations  of  the  Use  of  the  Graphs.  To  illustrate  the  use 

of  the  curves,  let  us  suppose  that  we  have  a  10,000  pound  (  m  ) 

maimed  interceptor  vehicle,  and  desire  to  intercept  a  6000  n.m. 

minimum-energy  warhead  trajectory  at  a  point  0.45  radians  (  wx  ) 

down— range  from  the  apogee.  We  also  select  a  launch  site  that  is  0.23 

radians  (  <£  )  up-range.  By  entering  Bug.  6-10  at  -0.23  on  the 

m 
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Launch  Site  Position  Relative  to  Warhead  Apogee  <pLS  (radians) 

fff 

Fig.  6-11 

The  Effect  of  <Pls„  on  Th^  for  Various  Intercept  Points  on 
6000  n.m.  High  Angle  Trajectory 
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t*  1  r* 


HSISHilfti 


180 


120  Ir-tH1"*! --r :: 

! Time  Limitation 


Thrust  Limitation 


.  "g"  Limitation 


Possible  Launch  CH'iSr.Zjr^ 
Sites  and  Intercept  :::f 
Points  in  This  Area  f7:--::TrT':~ 


—•4  — »2  ,  0  .2  .4  *6 

Launch  Site  Position  Relative  to  Warhead  Apogee  <pL$  (radians) 

m 

Fig.  6-lla 

The  Effect  of  Time,  Thrust,  and  ngn  Limitations  on  Location 
of  Possible  Launch  Sites  and  Intercept  Points 
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abscissa,  proceeding  vertically  upward  to  the  curve  labeled  0. 45 > 
and  then  horizontally  to  the  ordinate,  we  find  that  Th^  =  105-  We 
also  see  that  the  maximum  acceleration  is  less  than  10  g‘s  and  that 

the  excess  time  is  between  300  and  600  seconds.  By  using  the  same 

.  * 

procedure,  and  Fig.  6-13 »  we  find  that  T  =  410  seconds.  The  first 
stage  boost  thrust  required  is 


As  a  second  illustration  let  us  suppose  that  our  restrictions 
* 

are  -  600  seconds  and  v  *  10  g's,  and  we  want  to  intercept 

a  6000  n.m.  high-angle  warhead  trajectory.  Then  by  referring  to  Fig. 
6-11  we  see  that  we  are  restricted  to  operate  in  the  area  between  the 
600  second  dotted  line  and  a  curve  joining  the  stars.  This  area 
defines  the  possible  combinations  of  intercept  point  and  launch  site, 
and  the  corresponding  values  of  Th^  that  meet  the  restrictions. 
Furthermore,  suppose  we  have  a  5000  pound  interceptor  vehicle,  and  a 
first-stage  booster  rocket  with  1,000,000  pounds  of  thrust,  giving  us 
a  Th.  of  200.  We  are  then  further  restricted  to  launch  sites  and 
intercept  points  below  this  horizontal  line  on  Fig.  6-11.  Figure 
6-lla  shows  how  these  restrictions  can  limit  the  choice  of  intercept 
points  and  launch  sites. 


95 


GAW  Mech  61-2 


By  comparing  Pigs.  6-9,  10 ,  and  11  it  is  seen  that  as  we 

progress  from  the  low  angle  to  the  minimum  energy  to  the  high  angle 

warhead  trajectories,  time  and  g-limitations  become  less  restrictive 

hut  higher  values  of  Th^  are  required.  The  curves  of  Tx  (Pigs.  6-12, 

13,  and  14)  are  also  useful  for  determining'  the  possible  launch  sites 

and  intercept  points  for  a  given  time  restriction.  If,  for  example, 

we  detected  an  enemy  missile  on  a  6000  n.m.  minimum  energy  trajectory 

immediately  after  it  was  launched,  and  we  required  a  minimum  of  700 

seconds  to  launch  our  interceptor  system,  then  we  are  restricted  to 

use  values  of  4>  and  on  the  portions  of  the  curves  above 

~  m 

a  horizontal  line  at  700  seconds  (Pig.  6-13). 

Figures  6-15  through  6-20  show  basically  the  same  information 

for  the  4000  mile  warhead  trajectories  that  was  shown  in  Figs.  6-9 

through  6-14  for  the  6000  mile  trajectories.  The  dotted-line 

* 

isochrones  are  for  equal  to  300  and  450  seconds.  The  absence  of 

stars  on  some  of  the  curves  indicates  that  10  "g's"  was  never  exceeded 

for  these  cases.  These  curves  show  three  important  differences  for 

the  shorter  range  warhead  trajectory  :  maximum  acceleration  is  not 

as  restrictive;  lower  thrusts  are  required;  and  time  is  much  more 

critical.  Figures  6-18,  19  and  20  show  that  for  some  intercept 

* 

points  near  the  warhead  apogee,  T  becomes  negative.  This  means  that 
the  interceptor  flight  time  is  greater  than  the  warhead  flight  time, 
•and  that  the  interceptor  must  be  launched  before  the  warhead  for  these 
cases.  Since  this  is  quite  impractical  we  can  say  in  general  that 
for  shorter  range  enemy  missiles,  the  intercept  points  must  be  located 
further  down  range,  than  for  longer  range  trajectories.  It  is  also 
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The  Effect  of  <f>Lsmon  Tx  for  Various  Intercept  Points  on 
6000  n.m.  High  Angle  Trajectory 
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reasonable  to  assume  that  the  interception  of  an  8000  mile  missile 

would:  impose  more  severe  g-1 imitations;  require  higher  values 

of  first-stage  thrust;  and  increase  the  excess  time. 

Summary  of  Results 

1.  The  thrust  to  initial  weight  ratio  of  the  velocity-match 
rocket,  |3^  has  a  negligible  effect  on  overall  thrust  requirements. 

2.  0  has  an  appreciable  effect  on  maximum  accelerations  for 
trajectories  with  small  semi-range  angles. 

3-  Pu  has  a  negligible  effect  on  excess  time. 

4.  There  is  an  optimum  value  of  thrust  to  initial  weight  ratio  of 
the  boost  sub-rockets,  |3 ,  for  each  intercept  condition  which 
minimizes  overall  thrust  requirements. 

5-  For  values  of  p  close  to  the  optimum  value,  maximum 
accelerations  are  well  below  human  tolerance  limits. 

6.  A  small  increase  in  0  causes  an  appreciable  increase  in  excess 
time,  especially  for  small  semi-range  angle  trajectories. 

7«  Small  increases  in  both  structural  factor  and  specific  impulse 
greatly  decrease  overall  thrust  requirements. 

8.  Specific  impulse  ha3  a  much  larger  effect  on  thrust 
requirements  than  structural  factor  does. 

9.  An  increase  in  specific  impulse  causes  a  negligibly  small 
decrease  in  excess  time. 

10.  Structural  factor  has  no  effect  on  excess  time. 

11.  Excess  time  decreases  as  the  launch  site  and  intercept 
point  are  moved  up-range. 
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Launch  Site  Position  Relative  to  Warhead  Apogee  <pLSm  (radians) 
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The  Effect  of  9LSmon  T  for  Various  Intercept  Points  on 
4000  n.m.  Minimua  Energy  Trajectory 
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Launch  Site  Position  Relative  to  Warhead  Apogee  4>USm  (radians) 

Fig.  6-20 

The  Effect  of  Tx  for  Various  Intercept  Points  on 

4000  n.a.  High  Angle  Trajectory 
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12.  Thrust  and  weight  requirements  increase  as  the  launch  site 
is  moved  down-range  and  intercept  point  is  moved  up  range. 

13*  Thrust  and  weight  requirements  increase  as  the  warhead 
trajectory  altitude  increases. 

14.  Excess  time  increases  as  warhead  trajectory  altitude  increases. 

15 •  Higher  thrusts  and  weights  are  required  to  intercept  longer 
range  trajectories. 

16.  Excess  time  decreases  as  warhead  range  decreases. 

17.  Maximum  acceleration  increases  as  warhead  range  increases. 

18.  Time,  thrust,  and  "g"  limitations  determine  the  possible 
intercept  points  and  launch  sites  for  a  given  warhead  trajectory. 

19.  The  validity  of  the  assumptions  made  in  obtaining  the  numerical 
results  was  verified  for  a  typical  intercept  trajectory.  This 
makes  the  graphs  a  useful  tool  for  obtaining  approximate  thrust 
and  weight  requirements  for  any  size  of  interceptor  vehicle, 
for  the  intercept  conditions  considered. 

Conclusions 

The  trajectory-match  intercept  technique  appears  to  be 

theoi-etically  feasible: 

1.  If  first  stage  rocket  engines  of  sufficient  thrust  to  boost  the 
interceptor  vehicle  onto  the  intercept  trajectory  are  made 
available.  Minimizing  the  weight  of  the  interceptor  vehicle  is 
a  major  factor  in  satisfying  this  requirement. 

2.  If  sufficient  warning  time  is  available  from  ballistic  missile 
radar  or  satellite  reconnaissance  vehicles. 
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3.  If  low  reaction  times  are  possible  for  the  interceptor  system. 

One  major  requirement  to  meet  this  restriction  is  the  use  of 
solid  or  storable-liquid  propellants. 

4.  If  guidance  systems  are  developed  which  can  meet  the  requirements 
of  this  intercept  technique. 

Therefore  it  is  recommended  that  further  study  of  thi3  technique 
be  made  for  other  warhead  trajectories  and  for  different  combinations 
of  the  boost  and  velocity-match  stages.  This  study  could  also  be 
extended  to  the  more  realistic  situation  of  a  non-planar  intercept 
over  a  rotating  earth. 

In  addition,  if  the  interceptor  is  to  be  manned,  a  study  should 
be  made  of  the  methods  which  could  be  used  for  recovery  of  the  vehicle. 
Two  possibilities  are  a  controlled  skip-glide  re-entry  or  a  parachute 
recovery. 

Even  if  all  of  the  above  requirements  are  met  for  this  system, 
the  final  decision  as  to  whether  or  not  this  intercept  technique  is 
the  answer  to  the  ICBM  defense  problem  must  be  based  upon  an 
operational  analysis,  in  terms  of  cost,  technical  and  military 
manpower  requirements,  logistics,  and  many  other  factors.  Thus, 
the  trajectory-match,  ICBM  intercept  technique  appears  to  be 
theoretically  feasible,  but  a  more  comprehensive  study  is  necessary  to 
determine  if  it  is  militarily  practical. 
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Appendix  A 
Computer  Programs 

The  computer  programs  used  in  this  analysis  are  written  in  IBM 
1620  Fortran  Language  (Ref  5)  for  use  on  the  IBM  1620  digital 
computer.  The  instructions  for  the  use  of  these  programs  are  given 
in  this  appendix.  The  actual  Fortran  Language  Program  follows  each 
set  of  instructions.  The  paper  tapes  for  these  computer  programs 
are  on  file  with  the  Mechanics  Department  of  the  Institute  of 
Technology,  Wright-Patterson  Air  Force  Base,  Ohio. 

The  values  of  tq  and  ^  used  in  these  programs  are  (Ref  7s 430): 

(average)  =•  3437*747  nautical  miles 

p.  -  62,628.  (nautical  miles)//  (second)2 

Program:  CM-3  General  Trajectory  Evaluation 

This  program  gives  the  values  for  r,  v,  and  y  for  any  desired 
position,  ,  on  the  Ballistic  trajectory.  It  also  gives  the  time  of 
flight  between  apogee  and  ^  . 

Operating  Instructions. 

1.  V/hen  the  computer  prints  the  number  (l),  the  inputs  are  4>  q 
and  of  the  trajectory. 

2.  V/hen  the  computer  prints  the  number  (2),  the  inputs  are  the 
general  angle  ,  which  defines  the  point  to  be  investigated, 
and  L  ,  which  allows  the  computer  to  automatically  increase 
<p  hy  A  <p  after  each  computation.  This  gives  the  option  of 
determining  the  values  for  regular  intervals  on  the  trajectory. 
However,  if  only  one  point  is  desired,  make  A  greater  than 
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cfc  ,  since  some  value  of  A  must  be  inserted  into  the  computer 
memory.  When  <f>  becomes  greater  than  4>  ,  the  computer  sets  <fi 
equal  to  q,  calculates  the  values,  and  then  returns  to  the  start  of 
the  program  (ready  to  accept  new  values  of  <f3  Q  and  yq). 

3*  The  results  are  typed  in  the  following  order:  ,  T*,  r,  y,  and  v. 

4-  To  change  either  or  A  at  any  time,  turn  "Sense  Switch  1"  on. 
The  computer  will  then  print  the  number  (l),  and  new  values  can  be 
inserted  in  the  memory. 


^  ‘  y  t»  l  \'*Z. 


'  r  3  T  r,S2  ^ 


U-62628. 

RO-3437.747 
5  M-1 
PRINT, M 
ACCEPT ,PO,GA 
IF(PO)  1,2,1 

1  E-1.0/(COS(PO)+(SIN(PO))*(COS(GA))/(SIN(GA))) 
HS-U*R0*(1 .0-E*(C0S(P0))) 

H-SQR(HS) 

W-1 .0-(E**2.) 

G* (H**3 • )/ ( (U**2 • )*W) 

P-(H**2.)/U 
Z«(SQR(W) )/ ( 1 ,0-E) 

0*2.0/ (SQR(W)) 

6  hV2 
PRINT, H  ‘ 

ACCEPT, X,0ELTX 

8  S-(-E)*(SIN(X))/(1 .0-E*(C0S(X))) 

Y»Z*( S I N( X/2 . ) )/ ( C05( X/2 . ) ) 

T«G*(-S+Q*(ATN(Y))) 

R«P/(1 .0-E*(C0S(X))) 

A-ATN(S) 

V-H/(R*(C0S(A))) 

PRINT,X,T,R,A,V 
I F ( SENSE  SWITCH  1)  6,7 

7  IF(X-PO)  4,5.4 
4  X-X+DELTX 

I F(X-PO)  8,5,3 
3  X-PO 
GO  TO  8 
END 
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Program:  CM-11  Combined.  Boost-Phase  and  Intercept  Analysis 

This  program  contains  the  intercept  and  trajectory  geometry, 
rocket  equations  of  motion,  and  boost-phase  approximation  equations. 
These  equations  are  combined  as  outlined  in  Chapter  IV  to  give  the 
complete  solution  to  the  intercept  problem. 

Operating  Instructions. 

1.  When  the  computer  prints  the  number  (l),  the  inputs  are  n, 


2. 

3. 

4- 

5- 
6. 


I  and  (3  . 

8  U 

When  the  computer  prints  the  number  (2),  the  inputs  are  3  and  k. 

When  the  computer  prints  the  number  (3),  the  inputs  are 

and  v 
'wo 

When  the  computer  prints  the  number  (4),  the  input  is 


wx 


When  the  computer  prints  the  number  (5),  the  input  is  <£> 


mo 


Now  the  computer  solves  the  problem.  The  results  (for  k  ■*  3) 
will  be  typed  in  the  order  listed  below: 


moa 

*b 

u 


'moa 

kt. 


4>- 


4> 


mla 


u 


Y LX  A 

X 


Yo  ^1  ^2  Y  3 

Th  Th,  Th_  Th- 

u  1  2  3 

7.  The  computer  will  pause  after  each  line,  and  the  start  button 
must  be  pushed  to  continue  the  computation. 

8.  When  the  computation  is  complete,  the  operator  has  several 
options  as  to  which  part  of  the  program  the  computer  will  go 
to  for  the  next  set  of  intercept  conditions  to  be  investigated. 
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options  are 

controlled  by  the 

sense  switches 

as  follows; 

Change 

SW  1 

SW  2 

sw  3 

^  mo 

Off 

— 

— 

d> 

^  wx 

On 

Off 

— 

(3  or  k 

On 

On 

Off 

or  v 

wo  'wo 

On 

On 

On 

Program;  CM-12  Runge-Kutta  Integration 

This  program  is  designed  to  numerically  integrate  the  boost-phase 
rocket  equations  of  motion  for  a  three  stage  rocket  as  a  check  on  the 
results  of  the  boost-phase  approximations. 

Operating  Instructions. 


1.  when  the  computer  prints  the  number  (l),  the  inputs  are  I  ,  P, 

s 


and  t,  . 
b 

2.  When  the  computer  prints  the  number  (2),  the  inputs  are  Th^,  A^, 
Th2,  A2,  Thy  and  Ay 

3.  When  computer  prints  the  number  (3),  input  t^  (duration  of 
vertical  flight  segment),  a  (angle  between  thrust  vector  and 
velocity  vector  during  tip-over  phase),  and  T  (flight  path  angle 
where  tip-over  is  complete  and  a  is  set  equal  to  zero). 

4.  When  computer  prints  number  (6),  input  initial  conditions  for; 
t  (usually  zero)}  <fi  (usually  zero);  r  (usually  zero); 

Ls  (usually  -  ^moa);  r  (usually  rQ) ;  and  yQ  (usually  TT^). 

3.  Turn  Sense  Switch  1  on. 

o.  When  computer  prints  number  (4),  input  At,  the  incremental  time 
interval  of  the  integration. 
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7.  When  computer  prints  number  (5)>  turn  Sense  Switch  1  Off. 

Then  input  At  (again),  and  h,  (the  factor  which  determines 
the  print-out  interval).  The  results  are  printed  out  for 
time  intervals  equal  to  hA  t,  and  are  in  the  following  order: 
t,  <p  ,  r,  and  v. 
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Fortran  Language  Program  CM-11 


D*3. 

TD«=1  .OE-12 
TR«1 .0 
U-62628, 

RE«3^37.747 
P I  E-3.141 5927 
M»1 

PRINT, M 
ACCEPT, F.Z.BU 
21  M«2 
PRINT, H 
ACCEPT, B,Q 
20  fV3 
PRINT, H 

ACCEPT, POW.GOW 

SPOW*SIN(POW) 

CPOW-COS(POW) 

SGOW»SIN(GOW) 

CGOW«=COS(GOW) 

AA«SPOW*CGOW/SGOW 

EW»1 ,0/ (CPOW+AA) 

HWS=U*RE*( 1 .0-EW*CP0W) 

HW=SOR(HWS) 

AB-HWS/U 
19  M«4 
PRINT, M 
ACCEPT, PWI 
4  SXA«SIN(PWI ) 
CXA*COS(PW! ) 

AC«1 .0-EW*CXA 

AG«EW*SXA/AC 

RI-AB/AC 

gi=atn(ag) 

CGI«C0S(GI ) 

VW 1 *HW/ ( R I *CG I ) 

16 

PRINT, M 
ACCEPT, PO 

ZZ*((RI/RE)-1 ,0)/ AG 
CPO-COS(PO) 
ZA«S0R(ZZ*ZZ+1 .0) 
CAZ-CPO/ZA 
SAZ*SQR(1 .0-CAZ*CAZ) 
AZ— ATN(SAZ/CAZ) 
ZB=ATN(ZZ) 

P l*AZ-ZB 
SPI-SINfPI ) 

CPI»C0S(PI j 
E«AG/ ( AG*CP I +SP I ) 
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HS-U*RE*(1 .0-E*CP0) 

H-SQR(HS) 

VI-H/(RI*CGI ) 

DV 1-VW1 -V  I 
GOU/  (RE*RE ) 

EV-Z*GC 
DLU-DVI/EV 
DU-EXP ( DLU ) 

W0U-(F-1 ,)*DU/(F-DU) 

WU-WOU-1 . 

TBU*EV*(DU-1 . )/ (GC*BU*DU ) 
QQ-EV*EV*(DU-1 .-DLU)/ (GC*BU*DU) 
AJ«(PIE/2.)-G! 

SAJ-S I N( AJ ) 

CAJ-COS(AJ) 

TAK-SAJ/( (QQ/RI )-CAJ) 
AK-ATN(TAK) 

DP 1-P I E-AK-AJ 
SAK-SIN(AK) 

R I A-R I *SAJ / SAK 
G1A*AK-(P I  E/2. ) 

CG I A-COS (G I  A) 

HA-VI*RI A*CGI A 
HAS«HA*HA 
REC-HAS/'J 

ET-(VI*V!/2.)-(U/RI  A) 

EAS-1 ,+(2.*ET*REC/U) 

EA-SQR(EAS) 

AL-U*RE 

AM-U*RIA 

CPOA- ( AL-H AS ) / ( AL*E A ) 

CP  I A= ( AM-HAS )/ ( AM*EA) 

SPOA-SOR(1 ,-CPOA*CPOA) 

$PIA— S0R(1 .-CP  I A*CP I  A) 
POA-ATN(SPOA/CPOA) 

PIA=ATN(SPI A/CPiA) 

AFL-POA-PI A+DPI 
TGO«EA*SPOA/ ( 1 ,-EA*CP0A) 
GO-ATN(TGO) 

PB«0.9*P0A 
DA-0.1 *POA 
SPB-SIN(PB) 

CPB-COS(PB) 

12  AN-1 ,-EA*CPB 

TGB»EA*SPB/AN 

GB-ATN(TGB) 

CGB-COS (GB ) 

GBR-GB-DA 

SB-SIN(GBR) 

CB-COS(GBR) 

RB-REC/AN 

VB«HA/(RB*CGB) 

VX=VB*CB 
8  DD-D-1 . 


nl 


Mj.'.K  at -2 


DL*LOG(D) 

UK»EV*Q 
UU»EV*UK 
QW-Q-1  . 

DB*D*B 

YY«DD*CB/ ( DB*DL ) 

YZ-SQRO  ,-YY*YY) 
AP*ATN(YY/YZ) 

A=AP+GBR 

SA=S!N(A) 

CA=COS( A) 

F  D*UK*DL*CA-VX 

1  F ( FD*F D-TD )  1  ,1  ,2 

2  DFA«DD*CA+YY*(DD-DL)*SA/YZ 
DFD-=UK*DFA/(D*DD) 
D*=D-(FD/DFD) 

GO  TO  8 

I  AU*=2  ,* ( DD-DL)+QW*DD*DL 
AR«UU/(2„*GC*DB) 

XOAR*CA*Ai2 

YOAR*SA*Ai^-AR*Q*DD*DD/DB 
0A=ATN(XC/ ( YC+RE)) 
SDA-SIN(DA) 

RBA*=XC/SDA 
i F ( RB-RBA)  3,7,9 

3  RC-RBA-RB 
GO  TO  10 

9  RC-RB-RBA 
10  IF(RC-TR)  7,7,11 

II  RB*RBA 

CPB-(1 .-REC/RB)/EA 
SPB*SQR(1 ,-CPS*CPB) 

PB=ATN  ( SB/CPB ) 

GO  TO  12 

7  TB«UK*DD/(GC*DB) 

GOO* GO 
PBB=PB 
PIAA-PIA 
AFL  l=*AFL 

PRINT, POA , GOO , PBB , P I AA , AFL I 

PAUSE 

DBB=DB 

DBU=DU*BU 

ALF-A 

PRINT, TBU,TB,DBB,DBU,ALF 
PAUSE 

AW* ( F-1 . )*D/ ( F-D) 

WO B* AW* AW 
AS*1 .-(1 ./AW) 

AU*B/ AS 
THU=BU*W0U 
I  F  ( Q-2 . )  13J3.14 
13  W1 B=W0B*AS 
W2B=  (l.'0B-W1B)*AS 


ns 


•  'AW  Viftf'b  ^>]  —2 


WO»WOB*WOU 
W1«W1B*W0U  * 

W2«W2B*WOU 
PRINT,WU,W0,W1  ,W2 
PAUSE 

thi=wo*b 

TH2-AU*W2 
PRINT, THU ,  TH 1 , TH2 
GO  TO  15 

14  WOB«WOB*AW 
W1  B*WOB*AS 
W2B»(WOB-W1B)*AS 
W3B» (WOB-W1 B-W2B)*AS 

wo«wou*wob 

W1«WOU*W1B 

W2»WOU*W2B 

W3-WOU*W3B 

PRINT, WU, WO, W1 ,W2,W3 

PAUSE 

TH1»WO*B 

TH2-AU*W2 

TH3»AU*W3 

PRINT, THU, TH1 ,TH2,TH3 

15  pause 

I F (SENSE  SWITCH  1)  1 ? , 1 6 

17  I F ( SENSE  SWITCH  2)  18,19 

18  I F ( SENSE  SWITCH  3)  20,21 

END 


Pi  -i 


PI-3. 1415927 

EL=LOG(10.) 

U-62628. 

FM-6080.2011 

RE-3437.747 

DIMENSION  A(4)tE(4),C(4),TH(3).Q(5).X(5).Y(5),S(3) 
AL-O. 

GC=U/(RE*RE) 

M=1 

PRINT, M 
ACCEPT, Z,B,TB 
EV«Z*GC 
ZB-Z/B 
M=2 

PRINT, M 
DO  3  K-1,3 
3  ACCEPT, TH(K),S(K) 

17  NV=3 
PRINT, M 

ArrCOT  T\/  A  A  AO 

rixy  w  ui  i  •  i  t  ^  nr»  ^  ni 

DD-O. 

AZ*SQR(0.5) 

A ( 1 )»0.5 
A(2)*1 ,-AZ 
A(3)-1 .+AZ 
A(4)-1 ./ 6. 

E(1  )-2. 

E  ( 2  )*=1  . 

E(3)-1. 

E(4)*2. 

C(1 )-0.5 
C(2 )-A(2) 

C(3)-A(3) 

C(4)«0.5 
DO  2  1-1  .5 
2  0( I )— 0* 

PRINT, M 
DO  80  1-1  ,5 
80  ACCEPT, Y( I ) 

ACCEPT, GM 
L=11 
K-1 
N-1 
M-4 

PRINT,  M 
ACCEPT, H 
J=1 

37  V-SQR(Y(3)*Y(3)+Y(5)*Y(5)*Y(2)*Y(2)) 

GO  TO  (52, 53, 54), N 

53  0=0. 

TF=GM-Y(4 ) 

N=3 
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GO  TO  4 
54  TA-TF+Y(4) 

GO  TO  70 

52  VV-V*FM 

HT-Y(5)~RE 

HA-HT-86.67 

HB-HT-6.0855 

HC-HT- 13.487 

HD-HT-25.658 

HE-HT-28.783 

HF-HT-43.586 

I  F (HA)  18,18,19 

18  RH-1 ./ ( EXP( EL*(2 ,624+.108l8048*HT ) ) ) 
GO  TO  20 

19  RH-1 ./ (EXP( EL*( 1 2. +.01 142854*HA) ) ) 

20  IF(HB)  21  ,22.22 

21  AM-. 1 8362-. 0040095 *HT 
GO  TO  40 

22  IF(HC)  23,23,24 
23  AM-. 15922 
GO  TO  40 

24  IF(HD)  25,26,26 

25  AM-.1 5922+#00l8602*HC 
GO  TO  40 

26  I F (HE )  27  ,27,28 

27  AM-,18186 
GO  TO  40 

28  IF(HF)  29,30,30 

29  AM-. 181 86-. 0028798 *HE 
GO  TO  40 

30  AM-. 13923 
40  AN-V/AM 

I F (AN-. 6)  31,31,32 

31  CD=.25 
GO  TO  41 

32  I F (AN-1 .6)  33,34,34 

33  CD* AN-. 35 
GO  TO  41 

34  IF(AN-2.9)  35,36,36 

35  CD-I ,73-.3*AN 
GO  TO  41 

36  CD— ,86 

41  D-.5*CD*RH*VV*VV*S(K) 

4  TA-AL+GM 
70  ST-SIN(TA) 

CT-COS(TA) 

SG-SIN(GM) 

CG-COS(GM) 

X(1)-1. 

GO  TO  (73. 74, 75), K 

73  TS-Y(1) 

GO  TO  76 

74  TS-Y ( 1 )-TB 
GO  TO  76 
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75  TS*Y ( 1 )-2,*TB 

76  X(2)«EV*(CT-D*CG/(TH{K)))/(Y(5)*(ZB-TS))-2.*Y(3)*Y(2)/(Y(5)) 
X(3)«EV*(ST-D*SG/(TH(K)))/(ZB-TS)+Y(5)*Y(2)*Y(2)-U/(Y(5)*Y(5)) 
X(4)«Y(2) 

X(5M(3) 

DO  1  l«1  ,5 

WZ«A(J)*(X( I }-E(J)*Q( I )) 

Y( I )-Y( ! )+H*WZ 
1  Q( I )“Q( t )+3.*WZ-C( J)*X( f ) 

I  F ( Y( 2 ) )  71,72,71 

71  TG*Y(3 )/ ( Y (5 )*Y(2) ) 

GM=ATN(TG) 

72  J-J+1 
IF(J-4)  37,37,5 

5  J-1 

IF(Y(1)-TV)  6,7,7 

7  IF(GM-AP)  8,8,9 

9  AL*=AA 
GO  TO  10 

8  AL“0, 

10  I F ( Y ( 1  )-TB )  6,11  ,11 

II  I F ( Y ( 1 )-2 ,*TB ;  97,98,98 

97  K»2 

IF(L-22)  95,6,6 

95  L-22 
GO  TO  99 

98  ir(Y(1)-3.*TB)  96,93,93 

96  K-3 

I F ( L— 33 )  94,6,6 
94  L-33 
GO  TO  99 
93  K«4 
L-44 

99  PRINT, L 

6  I F ( SENSE  SWITCH  1 )  12,13 

12  IV5 
PRINT, M 
ACCEPT, HG.HH 
H«HG 

HP*Y ( 1 )+H*(HH-l . ) 

13  IF(Y(1)-HP)  14.15,15 

15  PRINT, Y(1),Y(4),Y(5),GM,V 
HP«Y(1 )+H*HH 

14  IF(K-3)  50,50,16 

50  IF(N-I)  49.49.37 
49  IF(D-DD)  38,39.39 
39  DD-D 

GC  to  37 

3P  ;F(D/(TH(K))-.0001)  51,37,37 

51  r*~s2 
GO  TO  37 

16  I F ( SENSE  SWITCH  3)  3,17 
END 
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^moa 

a 

0.50623 

rad 

d> 

r  b 

at 

-.41069 

rad 

^mla 

* 

0, 18201 

rad 

&LZA 

a 

0.68583 

rad 

\n 

» 

32.498 

sec 

Ktb 

= 

319.87 

sec 

V 

- 

4.1637 

"g" 

V 

u 

a 

4.9178 

’’g" 

X 

= 

0.73021 

rad 

V 

u 

a 

0.76462 

*0 

a 

72.974 

a 

51.872 

r= 

15.000 

a 

4.3376 

Th 

u 

a 

5.2939 

fti 

m 

109.46 

Th, 

8 

31.654 

Th3 

= 

9.1534 

By  inputting  the  corrected,  semi-range  angle  (<^>moa)  and  launch 
angle  (Tmoa)  for  the  intercept  trajectory,  a  print  out  was  obtained 
for  points  along  this  trajectory,  by  using  program  CM- 3. 

An  interceptor  vehicle  weight  of  10,000  pounds  (m^)  was  chosen, 
and  by  using  the  above  data  the  boost-stage  thrusts,  stage  weights* 
and  overall  system  weight  were  determined  to  be: 

T2  =  1,094,600.  lb 

T2  -  316,540.  lb 
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T3  -  91,534.  lb 

MQ  =  729,740.  lb 

m1  -  518,720.  lb 

IDg  =  150,000.  lb 

=  43,376.  lb 


The  burning  times  of  the  three  boost 


stages  are  equal,  and.  are 


_ b  =  106.62  sec 

3 

The  total  impulse  can  be  calculated  by  multiplying  the  thrust 
of  each  stage  by  the  burning  time  of  each  stage. 


Stage  1. 

5-84 

X 

104 

ton 

seconds 

Stage  2. 

1.68 

X 

104 

ton 

seconds 

Stage  3. 

0.49 

X 

104 

ton 

seconds 

Total 

8.01 

X 

104 

ton 

seconds 

The  equations  of  motion  used  are  (Ref  2:17) 

r  -  r<t>2  =  -jq-(Tstn  6  -h  L  cos  r-Dr  sin  <r)  - 
rep  -/-  2r<fi—  sos  0  —  L  sin  J  ~Dr  cos  ?r) 


where  r 

4> 

M 

T 


radius  from  center  of  earth 

angle  from  launch  site  to  radius  vector 

instantaneous  mass  of  interceptor  system 

thrust 
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L  -  lift  (neglected,  for  this  investigation) 

«  drag 

(a  «  earth's  gravitational  constant 
6  >=  angle  between  thrust  vector  and  local  horizontal 

Y  ■=  angle  between  velocity  vector  and  local  horizontal 

These  equations  are  broken  down  into  four,  first-order,  simultaneous 
differential  equations  for  the  computer  program  (Ref  9'.  110-120). 

Some  estimate  of  the  drag1  is  needed  for  these  equations,  since 
the  shape  and  size  of  the  interceptor  system  is  unknown.  A  curve 
of  vs  M  (mach  number),  (Ref  2:43) »  fo r  a  cylindrical  rocket  with  a 
40°  apex  nose  cone,  was  approximated  by  straight  line  segments  which 
resulted  in  the  following  estimate  for  C^. 

Cd  =  0.25  0  <  M  <  0.6 

Cd  =  0.25  +  (M  -  0.6)  0.6  <  M  <  1.6 

Cd  -  1.25  -  (M  -  1.6)  1.6  <  M  <  2.9 

C,  «  0.86  2.9  <  M 

a 

The  drag  is  given  by  the  equation 

Dr=irC„p  vs  A 

where  e  aerodynamic  drag,  lb 

«  coefficient  of  drag 
d 

p  =  atmospheric  density,  slug  per  ft^1 

v  *  velocity,  ft  per  sec 

2 

A  ®  cross  sectional  area,  ft 
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Now  since  the  coefficient  of  drag  is  a  function  of  mach  number, 
and  mach  number  is  a  function  of  the  speed  of  sound,  we  need  an 
estimate  for  the  variation  in  speed  of  sound  with  altitude.  This 
was  obtained  by  fitting'  straight  line  segments  to  a  curve  of  a 
versus  b,  where  b  is  the  height  above  the  earth  (Ref  1:24).  By 
doing  this  we  find  that 

a  -  0.18362  -  0.0040095  b  0  <b  <  6.0855 

a  =  0.15922  6.O855  <b  <  I3.487 


a 

=  0.15922  +  0.0018602  (b  -  13.487) 

13-487  <  b  < 

25.658 

a 

-  0.18186 

25.658  <  b  < 

28.783 

Si 

=  0.18186  -  0.0028798  (b  -  28.783) 

oa  ”70 ->  /  V  S 

|VJ  X  U  X 

A  1  c 

a 

=  0.13923 

43.586  <  b 

where  a 

is  the  speed  of  sound  in  nautical  miles 

per  second,  and  b  is 

the  height  above  the  earth  in  nautical  miles. 

In 

a  similar  manner  a  curve  of  log^  p  versus  b  (Ref  l:2l) 

is 

approximated  by 

log10  p  -  -  (2.624  +  0.1082  b), 

b  <  86.67 

log10  p  «  -(  12.0  +  0.01143  b), 

b  >  86.67 

The  remaining  factor  in  the  drag  equation  is  A,  the  cross 
sectional  area  of  each  stage.  We  must  obtain  estimates  for  all  three 
stages  since  burned-out  stages  are  jettisoned.  To  do  this,  first  we 
assume  that  the  propellant  density  is  0.06  lb  per  in  .  This  is  an 
average  of  densities  for  some  commonly  used  solid  propellants 
(Ref  10:  312,  313)»  With  this  density,  and  the  data  on  stage  weights, 
the  approximate  volume  of  each  stage  is  found  to  be: 
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Stage  1.  .  .  5 >000.  ft3 

Stage  2*  •  •  1  j445*  ft 

Stage  3»  .  .  420.  ft3 

With  the  volumes  known  we  can  now  assume  some  dimensions  for  the  stages 
and  estimate  the  cross  sectional  areas.  These  dimensions  are  as 
follows: 


Selected 

Approx. 

Approx. 

Area-ft^ 

Diam.  ft 

Length  ft 

Stage  1 

75 

10 

67 

Stage  2 

50 

8 

29 

Stage  3 

40 

7 

11 

The  velocity-match  rocket 

and  interceptor  vehicle  are, 

of  course, 

mounted  on  top  of  the  three  boost  stages  but  it  is  assumed  that  their 
cross  sectional  areas  are  equal  to  or  less  than  the  third  boost-stage 
area. 

The  approximate  equations  for  speed  of  sound,  coefficient  of 
drag,  and  atmospheric  density  are  built  into  computer  program  CM-12 
(Appendix  A),  along  with  the  equations  of  motion.  With  this  program, 
and  the  calculated  rocket  parameters,  a  trial  and  error  method  was 
used  whereby  different  values  of  the  launch  parameters  were  selected. 
The  values  of  r^,  v^,  and  were  then  compared  with  points  on  the 
desired  intercept  trajectory  printout,  to  determine  new  trial  input 
values.  After  several  trials,  the  following  results  were  obtained: 

Launch  parameter  inputs 

t  *  24  seconds 

v 
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a  •>  -0.0&  radian  (4*6°) 

T  *  0.8  radian  (45*8°) 


where  t 


is  the  duration  of  the  vertical  flight  segment 

is  the  angle  between  the  velocity  vector  and  the  thrust 

vector 

is  the  angle  where  tip-over  segment  ceases,  and  the 
zero-lift  (thrust  along  the  velocity  vector)  segment 
begins.  That  is,  when  y  reaches  a  value  of  0.8  rad, 
a  is  set  equal  to  zero. 


Computer  printouts 


CM-12  (actual) 

CM-3  (desired) 

Approx, 

Error 

rb 

3637.12  n.m 

3639.22  n.m. 

0.065C 

*b 

0.44583  rad 

0.43464  rad 

2.5* 

vb 

3.0322  n.m.  per  sec 

3.1342  n.m.  per  sec 

3.3£ 

The  approximate  error  in  overall  weight  (or  thrust)  due  to  the 
error  in  required  burnout  velocity  can  be  found  as  follows.  It  can 
be  shown  that  for  a  k-stage  rocket  in  field  free  space, 


r  y«- 

H,  r  n-D,e  *« 


Lfo-O,)'** 


where  is  the  overall  mass  (or  weight)  ratio  corresponding  to 


is  the  sub-rocket  effective  mass  ratio  required  to  reach 
burnout  velocity  v^, 


129 


GAW  Uecfc  61-2 


^  and  D-  correspond  to  burnout  velocity  v,_. 

^ 

If  we  use  the  preceding  computer  results,  and  the  above  equation,  with 

>  n  ■  10 


v. 


V, 


u 


1  -  -  2,7758 

{3  1*5 

3.1342  n.m.  per  second  (desired  burnout  velocity) 
3.0322  n.m.  per  second  (actual  burnout  velocity) 

!  g  .  (250X32.2)  ■  1.325  n.m.  per  second 


it  is  found  that 


h 

Tt 


0.8987 


actual  weight 
required  weight 


This  indicates  that  there  is  an  approximate  error  in  the  overall  Bystem 
weight  of  11$.  Since  the  first  stage  thrust  is  simply  0  times  the 
overall  weight ,  this  means  that  11$  more  thrust  is  required  to  reach 
the  desired  burnout  velocity,  or  approximately  1,220,000  lb  instead  of 
1,094,600  lb. 


It  is  felt  that  after  several  more  trials,  a  more  nearly  optimum 
flight  path  could  be  determined,  and  that  the  velocity  error  could  be 
reduced.  However,  these  results  indicate  that  the  assumptions  used 
in  the  mathematical  developsient  of  Chapter  IV  do  not  invalidate  the 
numerical  results  of  Chapter  VI.  In  any  case,  for  heavy  interceptor 
vehicles,  where  thrust  is  probably  much  greater  than  drag,  the  graphs 
in  Chapter  VI  can  he  used  to  determine  excellent  estimates  of  the 
thrusts  required  for  different  interoept  conditions. 
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